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ABSTRACT 

Multiyear civilian manned missions to explore the surface of Mars are 
thought by NASA to be possible early in the next century. Expeditions to Mars, 
as well as permanent bases, are envisioned to require enhanced piloted 
vehicles to conduct science and exploration activities. Piloted rovers, with 30 
kWe user net power (for drilling, sampling and sample analysis, onboard 
computer and computer instrumentation, vehicle thermal management, and 
astronaut life support systems) in addition to mobility, are being considered. 

The rover design, for this study, included a four car train type vehicle 
complete with a hybrid solar photovoltaic/regenerative fuel cell auxiliary power 
system (APS). This system was designed to power the primary control vehicle. 
The APS supplies life support power for four astronauts and a limited degree of 
mobility allowing the Primary Control Vehicle to limp back to either a permanent 
base or an accent vehicle. The results showed that the APS described above, 
with a mass of 667 kg, was sufficient to provide life support power and a top 
speed of five km/h for 6 hours per day. It was also seen that the factors that had 
the largest effect on the APS mass were the life support power, the number of 
astronauts, and the PV cell efficiency. From the results of the auxiliary PV mass 
study several points can be made: 

(1) The life support power requirements significantly affect the mass 
and size of the PV panels. For this reason, every effort should be 
made to design an energy efficient life support system and attempt to 
reduce astronaut requirements if possible, 
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(2) In addition to the life support power, the PV panel efficiency and the 
number of astronauts on board the rover also play a major role in 
determining the size and mass of the PV panels for the auxiliary power 
system. 

(3) Parameters that affect the mass of the auxiliary PV power system to a 
lesser extent, but are still important to the overall mass and size 
determination are the cruising speed of the PCV, and the surface 
density of the PV cells. 

(4) The parameter that has an insignificant effect on the auxiliary power 
system mass, but strongly impacts the volume of the fuel cells, is the 
reserve capacity of the RFC's. 

(5) In order to minimize the size and mass of the auxiliary power system, 
traversing should be done at night. 

The primary power requirement of 100 kWe was met using an SP-100 
type nuclear reactor coupled to a dynamic conversion system and a man-rated 
radiation shield. The masses of the shield, reactor, and the radiator were 
directly dependent on the reactor power, while this power was dictated by the 
electrical output and the conversion system parameters. The following four 
conversion systems were investigated for use in the PPS: Free Piston Stirling 
Engine (FPSE), He/Xe Closed Brayton Cycle (CBC), C0 2 Open Brayton Cycle 
(OBC), and Thermoelectric conversion. Of the systems studied, the FPSE and 
CBC systems yield the lowest overall PPS mass, with the Stirling system (15-16 
tonnes) giving slightly lower system mass than the Brayton system (18-19.5 
tonnes). Because the shield was such a large portion of the overall system 
mass, its size must be kept as low as possible. Therefore the system with the 
highest conversion efficiency has the lowest overall mass. 

Since the efficiencies calculated for both the Stirling and Brayton cycles 
were projected upper limits a parametric analysis was completed to see what 
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effect lower efficiencies would have on the total rover mass. The analysis show 
that while the masses of each system increased by reducing the efficiencies to 
29% for the Stirling and 25% for the Brayton, both systems still had a total mass 
that is significantly lower than the thermoelectric and the open Brayton systems. 

It was seen that major mass component for the PPS was the radiation 
shield. While it was possible to design an effective man-rated radiation shield 
(satisfying the dose rate limit of 0.3 Sv/y to the rover crew) for a manned Mars 
rover using a single layer each of tungsten and lithium hydride, the shield would 
not be mass efficient. A double layered W and LiH shield was more effective 
and significantly lighter. The thickness and positions of the tungsten layers, 
within the shield, were varied to determine the impact of secondary gammas 
produced in the W on the dose rate and total shield mass. 

The combined shield of LiH and W/depleted U was slightly heavier than a 
double layered W-LiH shield. The mass of the optimized man-rated W-LiH 
shield for the Mars rover vehicle powered by a nuclear reactor system increases 
with the thermal power of the reactor raised to the 0.3797 power, and varies 
from 8600 to 20580 kg over the power range of 100 to 1000 kWth, respectively. 
These shield masses were based on a conservative design correlation 
developed based on the results of 1-D neutronics analysis. Results show that 
shield mass estimates could be as much as 15 percent higher than those based 
on 2-D neutronics calculations. 
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1. INTRODUCTION 


1.1 Background 

For centuries the red planet has intrigued man. The nearest outward 
neighbor to the earth, Mars, has been the focus of much speculation as the 
home for alien life, although today it is known that the planet houses no 
intelligent life forms. The first successful probe to Mars was Mariner 4, which on 
July 14, 1965 flew to within 10,000 kilometers of Mars. In 1969 two more 
Mariner missions were sent to fly by Mars to photograph the surface. Mariner 9, 
the final Mariner mission, was sent in 1971 to photograph the two Martian 
moons, Phobos and Deimos, and perform a photographic survey of the Martian 
surface. This survey set the stage for the first two successful landing missions 
on Mars. On July 20, 1976 the Viking I lander touched down on the Chryse 
Plain on the surface of Mars. A month and a half later the Viking II lander also 
came to rest on the surface of our sister planet [Ezell and Ezell, 1984]. Although 
both the Viking missions were highly successful, NASA has not revisited Mars. 

The first step in the return to Mars may be an unmanned Mars rover 
Sample Return (MRSR) mission. A possible MRSR mission scenario would 
place two spacecraft in orbit about Mars and land a small (<1 000 kg) robotic 
rover on the surface of the planet. The rover's power would be supplied by one 
or two radioisotope thermoelectric generators [Shock, 1989]. This rover would 
collect up to 100 kg of soil, rocks, and atmospheric samples over a period of 1 1 
to 1 8 months. These samples would in turn be transferred to an ascent vehicle 
for return to Earth. In addition to collecting samples for return to earth, the rover 
would be capable of performing additional experiments over the lifetime of its 
extended mission. The Solar System Exploration Committee (SSEC) of the 
NASA advisory Council is strongly recommending that a Mars Sample Return 
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mission be undertaken before the end of the century [Bents, 1989]. This 
unmanned MRSR mission would then pave the way for manned exploration of 
our neighbor planet by plotting the terrain and studying, in depth, the 
environment which will be encountered. 

An evolutionary strategy is currently planned for manned missions to Mars. 
This strategy begins with an initial exploration of Mars and builds towards the 
establishment of an in situ propellant facility and ultimately a permanent 
manned base. Some of the mission design guidelines include an initial crew 
size of four growing to eight, reusability of selected vehicle elements, and 
aerobraking at both Mars and Earth. The program is to be carried out through a 
series of unmanned cargo flights and manned exploration missions [Andrews, 
1989]. 

1.2 Power System Options 

Multiyear civilian manned missions to explore the surface of Mars are 
thought by NASA to be possible early in the next century. Expeditions to Mars, 
as well as permanent bases, are envisioned to require enhanced piloted 
vehicles to conduct science and exploration. Piloted rovers, with extended 
range, have been identified as a viable means of achieving global access of 
Mars. For these missions, a rover vehicle with 30 kW e user net power is being 
considered. The operations covered by this power include: drilling, sampling 
and sample analysis, onboard computer and computer instrumentation, vehicle 
thermal management, and astronaut life support systems. In addition to the 30 
kW e user power, electric power will be needed to drive the rover across the 
Martian terrain. The only power system option that can provide this level of 
power at a reasonable mass penalty is a nuclear reactor heat source with a 
static or dynamic conversion system. Fore the auxiliary power system, which 
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will by used in the event of a primary power system failure there are several 
possible system options. 

1.2.1 Auxiliary Power System Options 

These options include photovoltaic (PV) collection systems, radioisotope 
thermoelectric generators (RTG’s), dynamic isotope power systems (DIPS), and 
high performance fuel cells. The following is a brief discussion of each of these 
systems. 

1 .2.1 .1 Photovoltaic Power Systems 

Photovoltaic (PV) power systems have been flown by NASA since the late 
1950's with the launch of the Vanguard satellite in 1958. Those first cells, while 
only delivering milliwatts of power (used to operate a tracking oscillator), 
operated for over six years demonstrating the reliability the PV systems [Ralph, 
1989 ]. The cells used for this mission were p on n silicon cells, which exhibited 
~5% efficiency at air mass zero (AMO). After 1962, NASA switched to n on p 
type silicon cells in order to improve the radiation resistance of the PV’s. During 
the 60's the efficiency of PV cells was increased to 12% by the introduction of 
gridded front contacts. This substantial increase in efficiency was further 
improved during the 1970's and 80's which saw cells improved to their current 
level of 15% AMO [Ralph, 1989]. Although efficiencies of 18.1% have been 
measured in the laboratory with simulated AMO conditions, it is predicted that an 
efficiency of 22% can be obtained with the use of techniques such as light 
trapping and surface passivation [Ralph, 1989]. The efficiency has increased 
over the years as has the resistance of the cells to radiation damage. 

Although Silicon has been the primary material used for the construction of 
photovoltaic cells in America's space program (due to the experience with the 
material, ease of manufacturing, and low cost of materials), several other 
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materials are now under consideration. The use of GaAs as Photovoltaic cell 
material also has certain inherent advantages. By using GaAs typical 
efficiencies for single junction cells of ~19% are currently being seen. The 
maximum calculated efficiency for GaAs cells is 27.5%. By going to dual 
junction cells of GaAs on Ge substrate, efficiencies of up to 22% have been 
obtained in the laboratory [Ralph, 1989]. For these cascade cells a maximum 
calculated efficiency of 35.7% is predicted [Landis et al., 1989]. While Gallium 
Arsinide cells have been used for many years in terrestrial applications, 
production for the space program has only recently begun. The primary 
disadvantage of these cells is their high cost. 

Another emerging technology in the area of space photovoltaics is the use 
of ultrathin films. Although there are several materials being researched, the 
primary candidates are amorphous silicon (a-Si) and copper indium diselenide 
(CulnSe 2 )- Thin film cells have several inherent advantages in that they have a 
high radiation tolerance, high specific power, can be formed into flexible 
blankets, and have a large manufacturing experience. However, thin films are 
handicapped by lower efficiencies, lack of spacecraft experience, and the fact 
that they are not currently produced on lightweight substrates [Landis et al., 
1989]. Although the spacecraft experience and light weight substraits are easily 
overcome by the increase in experience with the cells, the maximum efficiency 
is limited by the physical properties of the materials. Laboratory efficiencies of 
10% are now being produced for a-Si cells; however, by using multiple junction 
designs this could increase to greater than 15% in the future. Projected power- 
to-weight ratios for 10% cells should reach 350 W e /kg AMO in earth orbit [Ralph, 
1989]. 

The experience and information relating to silicon PV solar cells is quite 
extensive. However, while photovoltaic solar cells have been shown to be 
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extremely useful on orbital spacecraft, their effectiveness is directly dependent 
on the availability and intensity of the solar insolation. The solar insolation in 
turn is a function of the inverse squared distance the cell is from the sun. As 
one travels further from the sun the available insolation decreases. Since the 
orbit of Mars places it -1 .5 times further from the sun than the earth is, the 
average solar insolation is only 590 W/m 2 compared to the 1371 W/m 2 for an 
earth orbit [Appelbaum and Flood, 1989]. If this value is used with a projected 
efficiency of 30 %, the rover’s 30 kW e user power alone would require -340 m 2 
of panels. Although this area is already prohibitive it does not take into 
consideration the decrease in insolation due to atmospheric effects, nor include 
the energy losses due to storage for night consumption, and the power required 
by the rover for mobility. 

1.2. 1.2 Radioisotope Thermoelectric Generators (RTG’s) 

Radioisotope thermoelectric generators (RTG’s) are another possibility for 
an auxiliary power system for the manned rover. RTG’s are extremely reliable 
long term static power supplies which operate on the thermoelectric energy (TE) 
conversion principle. Thermal energy is produced internally by the natural 
decay of Plutonium-238 (ti/2 =87.7 yr) and is carried across an array of TE 
couples and rejected by radiative heat transfer through external fins [Kelly, 
1987]. The temperature difference between the heat source and the fins is the 
driving force for the TE conversion [Angrist, 1976]. RTG's have been 
demonstrated to be safe and reliable with millions of hours of flight experience. 
The first RTG (SNAP 3B) launched in 1961 was used to power the Navy’s 
Transit 4a and 4b navigation satellites. Although this small unit only produced 
4W e beginning of mission (BOM) power, it successfully operated for over six 
years. 
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The next generation RTG's were developed under the SNAP program and 
showed a substantial increase in power and specific power. The SNAP 1 9 and 
27 RTG’s were flown on the Viking and Apollo missions and produced a specific 
power of < 2.2 W e /kg. The SNAP RTG's employed lead-telluride (Pb-Te) 
thermoelectric couples for energy conversion [Hartman, 1988]. Pioneer 10 and 
1 1 missions are also each powered by four SNAP 19 RTG’s. Pioneer 10 has 
passed the mean orbit of Pluto and should have sufficient power to transmit 
data through 1990 (18 years after being launched) [Skrabeck, 1987]. Again, 
these units operated in excess of the predicted performance. The Multi- 
Hundred Watt (MHW) RTG was the next step in RTG design. These systems, 
which were used to supply power to DOD’s LES 8 and 9 spacecraft and 
NASA's Voyager 1 and 2, yielded a specific power of ~4.0 W e /kg (BOM). The 
MHW RTG used silicon germanium thermocouples (unicouples) for energy 
conversion and modular Pu-238 fuel sphere packs for the heat source 
[Hartman, 1988], 

The current generation of RTG uses the modular General Purpose Heat 
Source (GPHS) and unicouple thermoelectrics. The GPHS-RTG developed 
under sponsorship of the Department of Energy, has a specific power of 5.3 
W e /kg in ground tests for both NASA’s 4.2 year Galileo mission to Jupiter and 
the joint NASA/ESA 4.7 year Ulysses solar polar mission. The GPHS-RTG has 
a mass of 55.9 kg [Bennett et al., 1986]. 

Under current design is the Modular-RTG (MOD-RTG). This next 
generation RTG has the distinguishing feature of true modularity. By varying the 
number of modules assembled, power levels ranging from 20-342 W e can be 
obtained. The heat source for the MOD-RTG is again the Pu-238 GPHS; 
however, the TE converters are improved SiGe/GaP multicouples (each 
multicouple contains 40 thermoelectric legs). MOD-RTG’s are predicted to have 
specific power up to 7.7 W e /kg [Hartman, 1988]. 
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Although RTG's are extremely reliable, they have some inherent 
drawbacks. First, the net conversion efficiency of a static thermoelectric 
converter system is in the range of 4-8%. This implies that for the 30 kW e 
specified user net power of the rover, a total of 500 kWj would be needed 
(assuming an optimistic 6% efficiency), requiring ~882 kg Pu-238. In addition to 
being in short supply at present, the use of Pu-238 in a rover RTG power supply 
is highly unlikely due to economical and safety considerations. 

Another RTG drawback is that heat must be rejected during the journey to 
Mars, because the energy produced by the radioactive decay of the Pu-238 is 
continuous. Also, the power level is always decreasing (~1 .3% per year 
electrical power loss due to fuel decay) [Kelly, 1987]. For these reasons and the 
low specific power, RTGs are not a viable alternative for the rover's auxiliary 
power source. 

1 .2.1 .3 Dynamic Radioisotope Power System 

Another system that utilizes the GPHS isotope heat sources is the Dynamic 
Radioisotope Power System (DIPS). This type of system can be designed to 
produce power in the 1 to 15 kW e range with potential growth up to 20 kW e . 
Currently the DIPS program is being conducted under joint DOE/DOD 
sponsorship [Bennett, 1988]. In DIPS the radioisotope heat source is coupled 
to a dynamic energy conversion system which provides a higher conversion 
efficiency than TE's. Three thermodynamic cycles have been studied for use 
with DIPS; Closed Brayton Cycle, Organic Rankine Cycle (ORC), and Stirling 
cycle. Schematic diagrams for the CBC and ORC are shown in Figures 1 and 2 
[Bennett, 1988], respectively. Although system efficiencies up to 5 times that 
obtained using RTGs are possible, currently projected specific mass, ~150 
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Figure 1 . Schematic of Closed Brayton Cycle for DIPS [Bennett, 1988]. 



Figure 2. Schematic of Organic Rankine Cycle for DIPS [Bennett, 1988]. 
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kg/kW e (6.6 W e /kg), [Angelo and Buden, 1985; Pearson, 1988] of DIPS may 
limit their usefulness in a large, high power Mars rover. Also, DIPS must still 
contend with the same power degradation and heat rejection during transit 
issues as RTGs. 

1.2. 1.4 High Performance Fuel Cells 

High performance fuel cells such as Hydrogen-Oxygen primary fuel cells 
are a very efficient way to store energy. The major drawback with these 
systems is that they are an energy storage device rather than a production 
system. Thus, as the mission length and the total energy storage increases, the 
mass of the reactant storage becomes prohibitively large. Although 
Hydrogen/Oxygen fuel cells are not an acceptable primary energy source for 
the rover, they will be considered as a part of an auxiliary backup system in 
case of primary power system failure. 

1 .2.2 Nuclear Reactor Power Systems for Primary Power 

Nuclear reactor power systems exhibit many highly desirable 
characteristics as a primary power source for a manned Mars rover. First, the 
reactor core will produce virtually no heat until it has been brought up to power 
after safely landing on the surface of Mars. Second, the power output of the 
system can be controlled and adjusted to the level needed. Unlike PV systems, 
the power and sire is not dependent on the distance from the sun and the 
rotation of the planet. Current efforts in the area of space reactors are focused 
on the NASA/DOD/DOE SP-1 00 Program. This program is currently following a 
development plan that will provide the technical data base and the verified 
analytical design codes needed to design, fabricate and qualify space reactor 
power systems to meet future mission requirements. 
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The SP-100 nominal reactor is a fast spectrum lithium cooled reactor 
power system being developed under a joint DOE/DOD/NASA agreement. The 
goal of the SP-100 program is to develop a space nuclear reactor power system 
in the 100 - 1000 kWe range (GE 1988). The reactor core is coupled to 12 
primary coolant loops, referred to as the primary heat transfer system (PHTS), 
each equipped with a thermoelectric-electromagnetic (TEM) pump for 
circulating the Lithium coolant. The heat is transported from the reactor core by 
liquid lithium to the thermoelectric (TE) power conversion assemblies (PCA) 
where it is partially converted to electricity. The residual the heat is then 
transported by circulating lithium in the secondary loops to the heat pipe 
radiator where it is rejected into space. Each integrated power module (IPM) 
consists of a primary coolant loop and a secondary coolant loop having its own 
radiator panel, TEM pump, and PCA. Both loops are thermally coupled in the 
TEM pump and the PCA. To accommodate the change in lithium volume during 
operation and initial thaw upon deployment in orbit, each primary and 
secondary loop is equipped with a void accumulator/gas separator and a 
bellows type accumulator, respectively. The separator function of the 
accumulator in the primary loop is to remove helium gas produced in the reactor 
core by neutron-lithium reaction. 

The nominal thermal power of the reactor is approximately 2.4 MWt and 
the system electric output is 1 10 kWe resulting in an overall conversion 
efficiency of -4.5%. The SP-100 system is equipped with a W-LiH shadow 
shield that provides for a neutron fluence of 1013 and a gamma dose of 0.5 
Mrad for a 7-year operation time. The SP-100 fuel elements are highly 
enriched uranium/nitride with a Nb-.99%Zr-.01%C (PWC-11) cladding. The 
coolant is liquid lithium, while the structure of the reactor is composed primarily 
of PWC-1 1 . The core has an active height of 34-cm with a radius of 1 7-cm. 
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The control system for the SP-100 consists of BeO reflector shutter panels. 
When the shutters are opened more neutrons escape reducing keff of the core. 
In this work the beryllium reflector shutter panels are replaced by rotating 
BeO/B4C control drums (see Figure 3). This is done because, unlike the SP- 
100 reference mission, which employs a shadow shield that is not man-rated, 
the reactor for the vehicle will require (as a minimum) an oriented 4-n shield to 
insure adequate protection for the crew. Such a shield configuration limits the 
usefulness of a shutter control system. The current design of the SP-100 yields 
a specific mass of approximately 40 kg/kW e (25W e /kg). This design includes a 
shadow shield with a cone half angle of 17°. This will significantly increase the 
specific mass of the system. Even with the shield requirements the reactor 
system is felt to be the best alternative to give the rover sufficient power for 
global access. 

However, application of the SP-100 reactor power system to a manned 
Mars rover will require the resolution of a variety of technical issues in order to 
efficiently and effectively integrate the system into the rover vehicle. These 
issues include minimization of shield mass and optimization of shield 
configuration, optimal integration of the heat rejection system with the rover 
configuration, thermal management, protection from the harsh Martian 
environment, and contingency systems. The objective of thesis is to answer the 
questions that arise from these issues. While Chapters 2 and 3 focuses on the 
rover design and assessment of power requirements for traversing, Chapters 4, 
5, and 6 will report the results of parametric analysis to: (a) identify the suitable 
conversion system/heat rejection combination, (b) the radiation shielding 
design and integration of nuclear reactor with the rover, and (c) mass 
optimization of the rover. 
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1 .2.3 Objectives 

This research investigates the utilization of an SP-100 reactor power 
system to fulfil the primary power requirements for a manned Mars rover. The 
integration of a nuclear reactor power system into a manned Mars rover (see 
Figure 1 ) will require the resolution of a variety of technical issues. These 
issues include minimization of the man-rated shield mass and optimization of 
shield configuration, integration of the heat rejection radiator with the rover 
configuration, protection of the crew from nuclear radiation while having the 
flexibility to sample the Martian surface, and selection of an auxiliary power 
system for an emergency return scenario. In an attempt to address some of 
these issues the research focuses on assessing the auxiliary and traversing 
power requirements for the rover vehicle. A parametric analysis is performed to: 
(a) determine the sensitivity of the total mass of the rover vehicle to the number 
of astronauts on board, the desired cruising speed and the scenario for an 
emergency return in case the primary power system fails; (b) identify the type of 
energy conversion system for minimizing the mass of the nuclear power system, 
including the man rated shield; and (c) optimize the configuration and 
composition of the man-rated radiation shield to reduce the radiation dose to 
the crew. 
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2. ROVER LAYOUT AND DESIGN 


2.1 Design Considerations for a Manned Rover 

The integration of a nuclear reactor power system with a manned Mars 
rover requires a detailed description of the design and functional requirements 
of the rover. These requirements included: 

1 ) a 7-year lifetime with intermittent operations, 

2) unlimited range, 

3) life support for 4 astronauts, 

4) an SP-100 type reactor coupled to a static or dynamic 
conversion systems, 

5) a net user power of 30 kW e , and 

6) materials consistent with 7 year life in Martian environment. 

It should be noted that the SP-100 reactor containment and piping are 
fabricated from Nb-IZr, such material is not compatible with the C0 2 Martian 
environment. Therefore, this issue will have to be addressed through the use of 
a special protective coating or other suitable means to protect the structure from 
being attacked by C0 2 . The material compatibility however is outside the scope 
of this study, hence it will not be addressed further.These requirements form the 
starting point for a more detailed description of the rover in order to determine 
the mobility power requirements, and hence, determine the electric power 
output of the nuclear reactor system. 

Currently, there is little information on the design of a multi-purpose global 
access rover for Mars; however, extensive investigations on moon based rovers 
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have been completed [16,17]. Since many of the design parameters for a lunar 
rover are similar to those of a Martian vehicle many parallels may be drawn. 

The first step in understanding the characteristics of the rover vehicle is to 
define the design capabilities and requirements of the rover. However, the 
Martian scenario exhibits several unique qualities. Three major considerations 
for the Martian rover are: (1) The extreme distance from the Earth demands the 
lowest specific mass in order to minimize the launch and transportation cost; (2) 
The rover power system should contain the highest degree of reliability and 
redundancy with capabilities for an adequate auxiliary power system in case of 
a complete malfunction of the primary power system, (3) The rover must be 
compatible with a Martian atmosphere composed primarily of CO 2 . 

2.2 Rover Configuration 

There is a wide variety of possible configurations for a planetary surface 
rover, and many parameters which can affect their mobility power requirements. 
The soil/vehicle interaction is the most critical power variable and is a function 
of traction type. Different traction types including wheels, tracks, and walkers 
have been studied. In the area of wheels alone there are several options 
available including rigid, pneumatic tires, wire mesh tires, elliptical wheels, 
hemispherical (and cone) wheels, and hub-less wheels [Bekker, 1969]. 

Typical rover configurations could vary from a large single vehicle to an 
articulated crank mounted wheeled vehicle, to a multi-car overland train 
[Cataldo, 1989]. In order to minimize radiation shield mass radiation sensitive 
equipment and the crew must be located a sufficient distance from the reactor. 
This important requirement is most easily met with the multi-car overland train 
rover option, which allows the crew and radiation sensitive equipment to be 
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FIGURE 3: Layout for Train-Type Manned Mars Rover 



located in the forward cars, while the reactor and radiation insensitive 
equipment and supplies are located at the rear of the train. The proposed 
configuration for the manned Mars rover is shown in Figure 3. It consists of four 
units: the Primary Control Vehicle (PCV), the Experimental Unit (EU), the 
Supply Car (SC), and the Reactor Car (RC). The function and the component 
mass breakdown for each of the units are discussed separately in the following 
sections. 

2.3 Primary Control Vehicle and Auxiliary Power System 

The Primary Control Vehicle (PCV) is a double wall, pressurized, climate 
controlled car. The PCV houses the driving, navigation, and general control 
systems for the rover vehicle. The basic car design accommodates four 
astronauts. The PCV also includes the sleeping quarters, galley, a work station, 
life support systems and an auxiliary power system. The basic exterior 
dimensions of the vehicle are 3x3x10 meters, which provides an interior volume 
of approximately 80 m^. A generalized configuration for the PCV is shown in 
Figure 4. 

2.3.1 Auxiliary System Layout 

In the case of primary power system (PPS) failure the immediate concern 
is the survival of the rover crew. Due to the extreme distance from the earth it is 
essential that the rover be equipped with an emergency power system to 
maintain life support systems and a degree of rover mobility for several days, in 
the event of a PPS failure . For this purpose a PV/regenerative fuel cell power 
system is employed, see Figure 5 (a and b). Figure 5a illustrates the system 
configuration during solar energy collection. Here the PV’s not only supply 
energy for the PCV, but also to electrolyze H 2 O into H 2 and O 2 for the 
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FIGURE 4 : Layout for Primary Control Vehicle 








(b) 

Night Configuration 


Figure 5: Auxiliary Power System Configuration 
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regenerative fuel cells (RFC’s). In Figure 5b the night configuration is shown. 
Here, no energy is being collected by the PV’s and total rover power is supplied 
by the RFC’s. 

Since the primary mission concern In the event of a PPS failure is shifted 
from the experimental objectives to the survival of the crew, all nonessential 
equipment and supplies are to be jettisoned. The scenario calls for 
disconnecting the PCV from the rest of the train. The PCV then travels at a 
reduced speed back to a preestablished base or the mission assent vehicle. 
The PCV is equipped with up to three days worth of life support power stored in 
the Hydrogen/Oxygen RFC; the total mass of the fuel cell system including 
reactant is 720 kg. This Figure is based on a crew of 4 astronauts, a life support 
power of 1 kW e /astronaut, and a system specific energy of 400 W e h/kg [Eagle 
Engineering, 1988]. For this storage a total reactant mass of 107 kg is required 
(12 kg H 2 and 95 kg 0 2 ). For reactant storage tanks pressurized to ~2MPa 
[Angrist, 1976] a total storage volume of 9.5 m 3 is necessary. 

2.3.2 Mass Estimates of Auxiliary PV Panels 

In addition to the fuel cells, the auxiliary power system will include a set of 
Photovoltaic (PV) arrays in order to recharge the fuel cells on a daily basis and 
provide for both life support and mobility power. The mass of the PV array is a 
function of the total mass of the PCV (this includes the PV array itself), the 
driving scenario (daily range and cruising speed), and the efficiency of the RFC 
cycle. 

The PV panels mass can be given in a general form by: 


Mpv = Paux-SMpv. 


( 1 ) 


where P aux is the auxiliary power and SM pv is the specific mass (kg/kWe) of the 
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PV array. The auxiliary power required from the PV's can be divided into power 
required for life support, P| S , and power required for mobility, P m . While the life 
support power is constant for a given number of astronauts, the power 
requirements for mobility will be dependent on the terrain, emergency cruising 
speed, daily range, and the total mass of the PCV according to the following 
equation: 


P m = x(8.0x1 0‘ 5 kWh/kmkg)M tpcv 'V. (2) 

The value of 8.0x1 0* 5 kWh/kmkg for the soil traction coefficients is based 
on the average energy consumption of the Apollo Lunar rovers [Cataldo, 1989]. 
The traction factor modifier x is a correction factor to account for the difference in 
soil conditions and the difference in gravitational acceleration between the 
Moon and Mars and V is the PCV return velocity in km/h. Although the 
calculation seems straight forward at first, it is important to note that the PCV 
mass also includes not only the values shown in Table 1, but also the mass of 
the PV array. 

The mass for the the PV panels, which is a function of several parameter, 
can be expressed as: 

M pv - S y pv [n a Pis[tc+{td-tcVrifc] + x(8.0x1 0‘ 5 kW hr/km kg)(Mpcv+M pv ) V(tm/rif c )] (3a) 
tc 


and by rearranging: 

M n a pis[tc^td-tcVTifc] + x(80x1 0' 5 kWhr/kmkg)Mpcv V(tm/nfc) 
PV [tc/SMp V * x(8.0x10' 5 kW-hr/km-kg) V(tm% c )] 


(3b) 
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Table 1 . Primary Control Vehicle Mass Breakdown 
Based on Four Astronauts 
(From EEI Report 88-188) 


Structure and Pressure Vessel 
Inner Shell 

Mass 

(kg) 

490 

Outer Shell 

500 

Other Structures 

200 

Insulation 

130 

Galley 

70 a 

Personal Hygiene 

90 a 

Emergency Equipment 

30 

Man-Locks 

230 

EMUs 

680 a 

Avionics 

90 

Environment Control and Life Support 

200 a 

Workstation 

40 

Drive Stations 

80 

Sleep Quarters 

500 a 

Experiments and Payload 

500 

Crew 

360 a 

Active Thermal System 
Radiator 

160 

Pump 

20 

Heat Exchanger 

50 

Piping 

100 

Refrigerant 

300 

Wheels and Locomotion 

300 

Fuel Cells 

740 

Total 

5860 


a Values will vary depending on the number of astronaut 
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Equation (3b) gives the mass of the PV panels in terms of the collection time, tc, 
the mobility time during both the day and night (t m ), the number of astronauts 
(n a ), the cruising speed (V), the traction factor modifier (x), and the specific 
mass of the PV's (SM pv ). The total mass of the PCV then becomes the sum of 
the tabulated masses in Table 1 plus that calculated for the PVs by equation 
(3b). 

Determining the total PCV mass then becomes an issue of accurately 
estimating the PV parameters and identifying the sensitivity with which each of 
these parameters affect the PV mass. For the base case scenario the mass of 
the PCV and the number of astronauts have already been set at 5860 kg (Table 
1 ) and four, respectively. Since the gravitational acceleration on Mars is 
approximately twice that of the Moon a value of 2.0 is used for x- This value is 
comparable to the estimated values for the rolling resistance of the Martian soil 
of 0.35 and a rover drive efficiency of 50% [AeroVironment, 1989]. Although the 
value of 2.0 forx is 15% lower than that suggested by AeroVironment, it is still 
considered a reasonable value. 

The time parameters including PV collection (t c ) and rover mobility (t m ) are 
dependent on the driving scenario assumed; however, for the base case a 
collection time of 12.3 h (12.3 h is the number sun light hours at 0 degrees 
latitude during both the Martian Aphelion and Perihelion) and a mobility time of 
6 hours were used. A discussion on the selection of these base case values 
will be presented later in this report. Although the cruising speed for the 
emergency return, like the mobility and collection times, will be determined by 
the specific return scenario, a maximum speed of 5 km/h was assumed. An 
efficiency of 60% was used for the RFC system [Eagle Engineering, 1988]. The 
final parameter, the Specific Mass of the PVs, was calculated using an average 
surface solar insolation, G, on the planet’s surface of 250 W/m 2 , a PV efficiency, 


22 



r|p V , of 20%, and a panel surface density, SMA, of 2.3 kg/m 2 as: 

SM pv = SMA/(G-Ti pv /1 000 W/kW). (4) 

The base case values for the variable in equations 3b and 4 are listed in Table 
2 . 


Table 2. Base Case Auxiliary Power System Variables 


Photovoltaic panel efficiency 

npv 

= 20% 

Average solar insolation 

G 

= 250 W/m 2 

Specific mass area of PV panels 

SMA 

= 2.3 kg/m 2 

Fuel cell reserve life support energy storage (PCV) 

tap 

= 72 h 

Fuel cell reserve life support energy storage (EU) 

tae 

= 72 h 

Fuel cell reserve life support energy storage (SC) 

tas 

= 24 h 

Emergency return speed 

auxspd 

= 5 km/h 

Power required for life support 

Pis 

= ikW/astro. 

PV collection time 

tc 

= 12.3h 

Day Mobility 

tmd 

= 0 h 

Night Mobility 

tmn 

= 6 h 

Efficiency of the fuel cell cycle 

“Hfc 

= 0.6 

Specific mass of the fuel cells 

SM fc 

= 2.5 kg/kw h 


From the Base case parameters, PV panels having a mass of 667 kg are 
required to power the PCV during the emergency return. This corresponds to a 
PV panel area of 290 m 2 . If the panels were dimensioned such that they were 
1 0 m wide (the same as the length of the PCV), a 29 m length of PV panels 
would be needed. While this is not felt to be impossible to accomplish, a 
substantial increase in panel size would limit the auxiliary power system's 
usefulness to supply power for an emergency return trip, due to the large 
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storage volume and difficulty in deployment which would be associated with 
large panels. It is, therefore, important to evaluate the effects of reasonable 
variation in the base case parameters on the mass and area of the PV array. 

In order to minimize the PV mass it is important that the collection time be 
maximized and the mobility scenario be optimized for minimum power 
requirement. Therefore, it is assumed that energy collection takes place for 
12.3 h, representing all daylight hours, during which the PV panels are 
deployed and the PCV is not moving. Table 3 list the dark periods (in hours) at 
several Northern latitudes for both the Aphelion and Perihelion planetary 
positions [Kaplan, 1988]. Given 24.62 h as the length of a Mars Sidereal day 
and an equatorial location, the collection time can be determined by: 

t c = 12.3 hr-t^ (5) 

When determining the total mass of the PV panels it is important to consider the 
advantages and disadvantages of choosing a particular emergency return 
mobility scenario. The operating region for equatorial deployment is shown in 
Figure 6. The four curves presented in this Figure represent PV mass values for 
the base case parameters given in Table 2. The effect, upon the PV mass, of 
varying the collection time (or mobility during the day) and the mobility during 
the night, t mn , is presented. Starting from the base case of 6 hours of night 
mobility (@5 km/h this represents a daily travel of 30 km/day) and a 12.3 h 
collection time (0 h day mobility) a Photovoltaic mass of 667 kg is required. It is 
recognized that driving during the night can place additional constraints on the 
rover crew, but it reduces the mass of the PV panels, therefore, the merits of 
nighttime travel are discussed. 

Keeping the travel time constant at 6 h, a PV mass of 1500 kg is needed if 
traveling is limited to daylight hours. This increase from 667 kg to 1 500 kg 
(more than twice the base case mass) is due to the shortened collection time for 
the PV panels, compounded by the increased storage and consequent energy 
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loss during energy storage and recovery from the fuel cells. In addition to the 
mass increase, the size of the array becomes prohibitively large. For example, 
by traveling during the day instead of night, the 290 m 2 surface area is more 
than doubled to 652 m 2 . Therefore, in order to effectively reduce the size and 
mass of the PV, night travel during an emergency return is required. 


Table 3. Martian Seasonal Information 


Northern Hemisphere Summer (Aphelion) 


Latitude 

Mean Radius 

Ratio 

Dark Period 

(degrees) 

(km) 

(night/day) 

(hours) 

0 

3393.0 

1.00 

12.31 

9 

3351.2 

0.910 

11.73 

21 

3167.6 

0.793 

10.89 

33 

2845.6 

0.670 

9.88 

44 

2440.7 

0.538 

8.61 

55 

1946.1 

0.361 

6.54 

65 

1433.9 

0.000 

0.00 

Northern Hemisphere Winter (Perihelion) 



Latitude 

Mean Radius 

Ratio 

Dark Period 

(degrees) 

(km) 

(night/day) 

(hours) 

0 

3393.0 

1.00 

12.31 

9 

3351.2 

1.10 

12.90 

21 

3167.6 

1.262 

13.74 

33 

2845.6 

1.493 

14.75 

44 

2440.7 

1.858 

16.01 

55 

1946.1 

2.768 

18.09 

65 

1433.9 

undefined 

24.62 


When the collection time is held constant as in the base case, and the 
night mobility time is increased to 12.3 h (this assumes instantaneous 
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deployment of the PV’s), the PV mass increases to 900 kg. This value is still 
considerably less than the 1500 kg for only 6 hours of day travel. Had the 
additional 6.3 h been driven during the day instead, a PV mass of -2000 kg 
would be needed, requiring a total PV area of 870 m 2 , or 3 times that needed 
for the base case. 

The PV mass optimization surface, in terms of nighttime travel and daytime 
collection at 0 latitude (or daytime travel * 12.3h-t c ), is shown in Figure 7. It is 
apparent that the lowest PV mass requirement is that corresponding to 0 hrs of 
night travel and 12.3 hrs collection time (0 hrs day mobility). This scenario 
could be used in the case where a rescue mission is launched from a Martian 
base or orbiting space craft. This PV mass of 490 kg corresponds to the PV 
panels needed only to maintain life support. By traveling along the night 
mobility axis to 6 h, the base case in Table 2 is again reached. It is important to 
note that by increasing the night travel from 0 to 12.3 h, the PV mass is 
effectively doubled, while increasing the day travel from 0 to 8 h, increases the 
mass of the PV panels by a factor of 5. For this reason the night driving 
scenario is preferred. 

To illustrate the typical duty cycle that can be used, a log of daily activities 
is presented in Figure 8 for the base case scenario (see Table 2). This diagram 
shows the activities performed for a full Martian day, during an emergency 
return. It is assumed that the RFC's are fully charged (three days life support 
power) at sunset. One full hour each is devoted to folding and deploying of the 
PV panels. Although the actual time needed for these activities will depend on 
the amount of automation, an hour was determined to be more than sufficient for 
each activity, should the automatic deployment mechanisms fail. The schedule 
also allows for three meals and eight hours sleep each day. The 6 h driving 
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FIGURE 7: Auxiliary Power Systems Operating Surface 
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FIGURE 8: Base Case Mission Scenario for Emergency Return @ 5km/hr 



scenario allows for a total traveling distance of 30 km/day. It is important to 
reemphasize that in the scenario delineated in Figure 6, the entire 12.3 h of 
daylight is used for solar energy collection. The daylight hours during which the 
crew is awake, are designated for sample collection and scientific duties. The 
power utilization curve shown above the duty cycle represents the average 
power requirements based on life support power of 1 kW e /astro and a 
maximum cruising speed of 5 km/h. In the case where daily range must be 
increased, the night driving time can be increased from 6 to 10.3 h. For this 
scenario the total PV required mass is 820 kg. This 20% increase in mass and 
area of the PV panels over the base case mass can increase the daily travel 
distance by 70% (from 30 km/h to 51 km/day). For a 1000 km return trip only 25 
days will be required as compared with the 33 days for the base case. A 
possible duty cycle for the 10.3 h daily driving time is shown in Figure 9. Here 
the time allocated to scientific duties is decreased from 6 hrs to 2.7 hrs in order 
to free up time for travel. In Figure 7, the average power utilization of the longer 
driving time scenario is shown above the duty cycle. 

Figure 1 0 is a plot of PV mass contours as a function of both the collection 
time and the nighttime travel time. Similar to the surface in Figure 7, this plot 
allows for all the night/day driving combinations to be examined for given a 
mass requirement for the PV panels. In Figure 10, the minimum night driving 
scenario, maximum night driving scenario, and the base case scenario masses 
are all labeled with their respective locations on the contour. The lowest 
PVmass of 490 kg, represents the mass of PV panels needed to only supply 
sufficient energy to run the life support systems. 
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NIGHT TRAVEL (hr) 

FIGURE 10: PV Mass Contours 






2.4 Parametric Analysis of Design Parameters on PV Mass 

Because of the conceptual nature of the rover design the exact values for 
the parameters to determine the auxiliary power system mass are unknown. 
Even though care has been taken to use values judged to be reasonable, given 
current technology and expected advances, a degree of uncertainty in the 
results is still introduced. For this reason a parametric analysis was performed 
to examine the sensitivity of the mass of the auxiliary power system to various 
operation and design parameters in equation 3b. These parameters include 
cruising speed, number of astronauts, life support power per astronaut, average 
PV specific power, efficiency and specific mass of the PV panels, and the 
amount of fuel cell reactants stored. The results of the parametric analysis are 
presented and discussed in the following subsections. 

2.4.1 Effect of Cruising Speed 

The easiest variable to control in equation 3b is the Primary Control 
Vehicle's cruising speed. While decreasing the cruising speed reduces the 
mobility power requirements for a given traveling distance, it increases the 
duration of the return trip. Therefore, the appropriate cruising speed should be 
selected based on considerations of power requirements and traveling range. 

In Figure 1 1 , the required mass of the PV panels is plotted against the the PCV 
cruising speed with life support (kW e /astro) as a parameter. The points 
corresponding to 0 km/h cruising speed show the PV mass required to only 
maintain life support. The effect of cruising speed on PV mass is seen to be a 
relatively weak function. 
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2.4.2 Effect of Life Support Power Requirements 

As Figure 1 1 shows, maintaining a constant cruising speed of 5 km/h, and 
doubling the life support power to 2 kW e , increases the PV mass from 667 kg to 
just under 1200 kg. The strong dependence of the PV mass on the life support 
power is due to two considerations. First, unlike the cruising speed the life 
support power is continuous throughout the day. Second, for a crew of four 
astronauts doubling the life support power per astronaut increases the total life 
support power requirement by a factor of 4. Further increase of the life support 
power to 3 kW e /astro will bring the PV mass requirement to nearly 1 700 kg for 
the base case. As demonstrated in Figure 1 1 , an auxiliary power system that is 
designed for 1 kW e /astro at a cruising speed of 10 km/h would not be able to 
deliver sufficient power to maintain the life support at 2 kW e /astro, even without 
mobility. These results suggest that great care should be taken to minimize the 
life support power requirements. 

2.4.3 Effect of PV Specific Power 

In Figure 12 the PV mass is plotted against the PCV cruising speed for 
different values of the PV specific power . The specific power values in Figure 
12 are simply the average solar insolation (G = 250 W/m 2 ) multiplied by the 
efficiency of the array (from 8 to 40%). As Figure 12 shows, lowering the 
specific power of the PV panels increases the rate at which the PV mass 
increases with cruising speed. This rate increase is caused by the fact that PV 
panels with lower specific power require a larger collection area of the PV 
panels, resulting in a larger mass. In turn, this increased mass requires more 
power to move. In Figure 12, beginning at the base case, increasing the 
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CRUISING SPEED (in km/hr) 

FIGURE 1 2: Effect of PV Specific Power on PV Mass 




specific power to 100 W/m 2 brings the PV mass down to 300 kg. Conversely, 
decreasing the PV specific power to 20 W/m 2 raises the PV mass to as much as 
1750 kg. The sensitivity of the PV mass to its specific power is much greater 
than that of the cruising speed at high specific power values (>50 W/m 2 ); 
however, as the specific power decreases the sensitivity of the PV mass to the 
cruising speed increases drastically. These results suggest that the savings in 
the PV mass by lowering the cruising speed becomes less significant as the PV 
panel becomes more efficient. Therefore, future advances in PV technology 
should serve to shorten the travel time on Mars for the same PV panel mass. 
Another factor that could effect the specific power of the PV panels are the 
frequent dust storms that occur on Mars. These storms can significantly reduce 
the amount of sunlight received by the PV array (thus effectively reducing G). 
This reduction in received light can easily reduce the specific power by a factor 
of 2 for several days. The loss of power from the PV arrays could, however, be 
compensated for by using the reserve in the fuel cell storage system. 

2.4.4 Effect of PV Panel Efficiency 

The Specific Power of the PV panel, SPpv, is related to the PV cell 
efficiency, r|p V , as: 

SPpy = G'T|p Vi (6) 

where G is the average solar insolation available at the surface of the planet. 

An expected range of 20 to 30% for cell efficiency is shown in Figure 1 3. It 
should be noted that the base case value of 20% was chosen because it 
represents the projected efficiency for amorphous silicon (a-Si) cells which are 
light, flexible, and easily stored. If the PV efficiency increases to 25%, the same 
PV panel mass needed for the base case cruising speed of 5 km/h could power 
the PCV at a cruising speed over 10 km/h. Moreover, increasing the efficiency 
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CRUISING SPEED (In km/hr) 

FIGURE 13: Effect of PV Panel Efficiency on PV Mass 







to 30% would increase the PV effectiveness to the point where a 10 km/h speed 
for 6 hours of mobility would require less PV mass than that required for life 
support alone at 20% efficiency. Therefore, future advances in PV technology 
are expected to have significant impact on reducing the PV mass. 

2.4.5 Effect of PV Surface Density 

Another parameter that is completely independent of the efficiency is the 
PV surface density, or the amount of mass per square meter of PV panels. The 
value used for this parameter in the base case is 2.3 kg/m 2 . This value is based 
on current silicon cell technology for a sun tracking system, with projected cell 
efficiency of 20%. With the a-Si cells the specific mass is expected to decline to 
as low as 1.0 kg/m 2 , resulting in a 200% reduction in the mass of the PV panels 
for the base case. While the mass of the PV cells decreases as the surface 
density goes down, the total PV area required does not decrease significantly. 
The PV panel area, however, decreases as the PV cell efficiency increases. 
Therefore, the trade off is whether to use a larger, less efficient, yet lighter a-Si 
panels, or smaller, more efficient, but heavier tracking GeAs panels. In either 
case, the values presented for the base case should be conservative. 

Figure 1 4 demonstrates the effect of a surface density for 1 .0, 1 .65, and 2.3 
kg/m 2 . Note that by using a surface density 1 .0 kg/m 2 , the total PV mass 
decreases from the base case (SD * 2.3 kg/m 2 ) 667kg down to as low as 
280kg, but the array area only decreases from 290 m 2 to 280 m 2 . This is 
because the mass of the PV panels represents a small fraction of the total PCV 
mass. 

2.4.6 Effects of the Number of Astronauts and Reserve Energy Storage 
The final two parameters that affect the mass of the PV system are the 

number of astronauts on board and the amount of reserve energy storage of the 
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regenerative fuel cells. While both affect the PV mass by increasing the PCV 
mass and mobility power requirements, the number of astronauts has a 
profound effect on the life support power requirements, which is such a 
significant portion of the total auxiliary power needs. In Figure 1 5 the mass of 
the auxiliary PV panels is plotted versus the PCV cruising speed for 2, 4, and 6 
astronauts with the base case life support power of 1 kW e /astronaut. As Figure 
15 shows, the change in slope with the number of astronauts is hardly 
noticeable, yet the offset between the three curves due to the difference in life 
support power requirements is significant. These results again emphasize the 
importance of employing an efficient life support system. As demonstrated in 
Figure 1 5, a PV system that provides for four astronauts at a cruising speed of 
10 km/h, is sufficient to sustain six astronauts at a reduced speed of ~3 km/h. 

Increasing the amount of reactants for the reserve fuel cells does not 
significantly change the mass of the PV panels, since the increase in the 
reactant mass only affects the mobility power requirement. Also, because the 
reactant mass is a small fraction of the total PCV mass, increasing the fuel cell 
capacity from 24 to 168 h (1 day to 1 week) has only a small affect on the PV 
panel mass (see Figure 16). Since the fuel cells on board the PCV are 
regenerative types and the reactants are not stored as cryogenic liquids, but as 
compressed gases, it is not their weight, but their storage volume that becomes 
prohibitive. For example, a one day supply of reactants requires a total storage 
volume of 3.1 m 3 (1 .05 m 3 for 0 2 and 2.05 m 3 for H 2 ), while for 72 h, the 
storage volume will be 9.5 m 3 (3.2 m 3 for 0 2 and 6.3 m 3 for H 2 ). For one week 
storage, a reactant volume of 22.1 m 3 will be needed (7.4 m 3 for 0 2 and 14.7 
m 3 for H 2 ). 
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CRUISING SPEED (In km/hr) 

FIGURE 16: Effect of Fuel Cell Reactant Storage Capacity on PV Mass 






To summarize the results of the auxiliary PV mass study several points can 
be made: 

(1) The life support power requirements significantly affect the mass and 
size of the PV panels. For this reason, every effort should be made to 
design an energy efficient life support system and attempt to reduce 
astronaut requirements if possible. 

(2) In addition to the life support power, the PV panel efficiency and the 
number of astronauts on board the rover also play a major role in 
determining the size and mass of the PV panels for the auxiliary power 
system. 

(3) Parameters that affect the mass of the auxiliary PV power system to a 
lesser extent, but are still important to the overall mass and size 
determination are the cruising speed of the PCV, and the surface 
density of the PV cells. 

(4) The parameter that has an insignificant effect on the auxiliary power 
system mass, but strongly impacts the volume of the fuel cells, is the 
reserve capacity of the RFC's. 

(5) In order to minimize the size and mass of the auxiliary power system, 
traversing should be done at night; the daytime travel should be 
limited to as short a time as possible. 

(6) The PV mass calculated for the base case scenario appears to be 
reasonable. Therefore, the mass of 667 kg was added to the PCV 
mass of 5860 kg from Table 1 , giving a total PCV mass of 6527 kg. 

2.5 Experimental Unit 

The Experimental Unit (EU) is similar to the Primary Control Vehicle in that 
it is a double walled pressurized car capable of maintaining a crew of four 
astronauts in a climate controlled environment (see Figure 3). This car provides 
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facilities for sample collection and laboratory analysis as well as additional 
living space for extended missions. The EU also houses the control systems for 
the primary power system (PPS). The components and equipment pertaining to 
the EU are listed along with their respective mass estimates in Table 4. Again, 
these masses are based on the Eagle Engineering Inc. Report on Lunar 
Surface Transportation Systems, but have been modified to suit the manned 
Mars rover mission. 

Since the EU is not equipped with a PV/RFC auxiliary power system. In 
the event of a PPS failure, the EU would be detached from the PCV and 
abandoned. In the event of a minor, short term, potentially repairable 
malfunction in the PPS, the EU is equipped with three days worth of life support 
power in the form of primary fuel cell (PFC). The reactants for these cells are 
stored in liquid form in order to decrease the required volume. These reactants 
are included to give the rover crew additional time during which they could 
initiate repairs, before they would be forced to begin an emergency return. 

2.6 Supply and Storage Car 

The Supply and Storage Car (SC) is much simpler than either the PCV or 
the EU. The SC is an open trailer which carries all the radiation hardened 
equipment and tools (see Figure 3). At the rear of the car, a thermal shield is 
installed to protect the contents of the car from the thermal radiation rejected by 
the waste heat radiator surface mounted on the reactor car. The SC also 
contains an additional 24 h worth of reactants for the PFC in order to 
supplement that in the EU. This yields a total of one weeks worth of reserve 
energy between the PCV, EU, and SC. The mass breakdown for the SC is 
given in Table 5. 
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Table 4. Experimental Unit Mass Breakdown 
(EE I Report 88-188) 


Structure and Pressure Vessel 

Mass 

(kg) 

Inner Shell 

490 

Outer Shell 

500 

Other Structures 

200 

Insulation 

130 

Galley 

70 

Personal Hygiene 

90 

Emergency Equipment 

30 

Man-Locks 

230 

Work Station 

40 

EMU’s 

340 

Reactor Control Systems 

100 

Avionics 

90 

Environment Control and Life Support 

200 

Showers 

80 

Experiments and Payload 
Active Thermal System 

900 

Radiator 

160 

Pump 

20 

Heat Exchanger 

50 

Piping 

100 

Refrigerant 

300 

Wheels and Locomotion 

300 

Fuel Cells (Reactants^ 

jm. 

Total 

5500 
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Table 5. Storage Car Mass Breakdown (EEI Report 88-1 88) 



MASS 


kg 

Thermal Shield 

200 

Wheels and Locomotion 

300 

Other Structure 

500 

Experimental Equipment 

887 

Fuel Cells (Reactants) 

240 

Additional SupDlies and Eouio. 

1000 

Total 

3127 


2.7 Reactor Car 

The reactor car houses the reactor systems, energy conversion systems, 
waste heat rejection equipment, and shielding for the protection of the crew. 

The mass of the Reactor Car (RC) cannot be determined by simply tabulating 
the components and summing their masses because mobility power depends 
upon the mass of the rover (including the RC), and the mass of the reactor 
system is a function of the total power needed. Therefore, an iterative approach 
similar to that taken for the PV auxiliary system, is needed to determine the total 
power required. 

A model was developed to determine the total electric power requirements 
of the PPS as a function of: mass of rover Utility Cars (RUC, including the PCV, 
EU, SC), cruising speed, number of astronauts, user net power, mobility traction 
parameter, and the specific mass of the PPS. The total electric power delivered 
by the reactor system is given by the general form: 

P, = PA + x(8.0x1 0* 5 kW*h/km*kg)M(-V. (7) 
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The mass of the reactor power system, which is proportional to its power output, 
can be given as: 

M r = SM r P t . (8a) 

The reactor power system mass includes the mass of the following 
subsystems: reactor, primary heat transport, shield, energy conversion, 
secondary heat transport and power management. In addition, the frame and 
drive system mass of the Reactor Car were also estimated at 1 0% of the total 
load mass. Hence equation (8a) becomes: 

M r = 1.1(SM r P t ). (8b) 

The total rover mass (Mj) is the sum of the masses of the individual rover cars: 

M t = Mpcv+M^+Mgc+1 .1 (SM r P t ). (9) 

Substituting equation (9) into equation (7) and rearranging results in the 
following general form for the total power requirement as: 

p x(8.0x10~ s kW h/km kg)[Mpcv»-MEu-»-Msc] V 
1 -x(8.8x10‘ 5 kW h/km kg) SM r V 

It can be seen from equation (10) that the total power requirement for the 
manned Mars rover depends on five primary variables: P a , %, SM r , V, and the 
masses of the Utility Vehicles (M pcv , M eu , M sc ). The next chapter presents the 
results of parametric analysis investigating the effects of each of these variables 
on the total rover power requirements . 
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3. PARAMETRIC ANALYSIS OF TOTAL MASS AND POWER 
REQUIREMENTS FOR MANNED MARS ROVER 


In this analysis, the values for the number of astronauts, the cruising 
speed, the surface traction parameter, power system specific mass, and the 
utility car mass were varied so that their effects on the total rover mass and the 
net electric power needed could be evaluated. The parameters for the base 
case are listed in Table 6. 

Table 6. Base Case Rover Variables 


Additional Power above mobility requirements 

^a 

= 30 kW e 

Mobility traction parameter 

X 

= 0.16Wh/kmkg 

Drive Efficiency 

^dr 

= 50% 

Specific mass of the Power system 

SM r 

« 100 kg/KW e 

Mass of the Primary Control Vehicle (w/o PV's) 

Mpcv 

= 5860 kg 

Mass of the Experimental Unit 

^eu 

= 5400 kg 

Mass of the Supply Car 

M sc 

= 3100 kg 

Minimum rover speed 

minspd 

= 0 km/h 

Maximum rover speed 

maxspd 

= 30 km/h 

Number of astronauts 

N a 

= 4 


3.1 Effect of Number of Astronauts 

Variation of the number of astronauts has a minimal effect on the total rover 
mass. Increasing the number of astronauts raises the mass of the PCV in Table 
1 , and also increases the auxiliary power system mass and mobility power 
requirements. This mass increase, however, is a small fraction of the total rover 
mass and does not significantly affect the mobility power. This means the 
power system mass is relatively insensitive to the number of astronauts. Figure 
17 shows that for a cruising speed of 20 km/h, increasing the number of 
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MAXIMUM CRUISING SPEED (km/hr) 
FIGURE 17: Effect of the Number of Astronauts on Total Rover Mass 




astronauts from 4 to 6, increases the total rover mass by only 8% (from 25,000 
kg to 27,000 kg). Most of this mass increase (-1600 kg) reflects the weight of 
the additional two astronauts, their gear, and the increase in the auxiliary power 
system mass to provide for the additional astronauts. The remaining mass (less 
than 400 kg) represents the incremental increase in the mass of the primary 
power system for meeting the increased mobility and life support power needs. 

3.2 Effect of the Cruising Speed and Soil Traction Parameter 

The maximum cruising speed of the rover could vary significantly with the 
terrain encountered; however, 30 km/h was selected as an upper limit. Since 
the gravitational acceleration on Mars is almost twice that on the moon, the 
weight of the rover vehicle force will increase proportionally and the vehicle will 
exhibit a higher rolling resistance. The wheel design and soil consistency may 
also vary, which will affect the mobility power requirement. It is, therefore, also 
necessary to study the effects of varying the soil traction parameter. 

The results delineated in Figure 17 show the strong dependence of the 
total rover mass on the cruising speed. The mass corresponding to a zero 
speed represents the mass of the rover vehicles, including that of the PPS, 
needed to meet the User power requirement (excluding mobility). As this figure 
shows, increasing the cruising speed increases the mobility power 
requirements, resulting in higher masses for the PPS and the rover vehicle. As 
Figure 17 shows, for the base case, the total mass of the rover is 26 metric tons 
at a cruising speed is 19 km/h. As the cruising speed increases to 30 km/h, the 
total rover mass for the base case of four astronauts increases to 32 metric tons 
(-12% increase), and the power requirements for the PPS increase to 160 kW e 
(a 60% increase) 

The effects on mobility power of varying the maximum cruising speed and 
soil traction parameter are illustrated in Figure 18. Here the mobility power 
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requirement is plotted against the maximum cruising speed, from 0 to 30 km/h, 
for three different surface traction values (8.0x1 0* 5 , 1.6x10* 4 , and 2.4x1 O' 4 
kW h/kg km). Starting with the base case of 1 .6x1 0* 4 kW h/kg km (x = 2.0), and 
a cruising speed of 19 km/h, the required power level is 70 kW e for mobility (P t 
= 100kW e ). Increasing the cruising speed to 30 km/h will require increasing 
mobility power 120% to 153 kW e . The effect of variation in the surface traction 
parameter becomes highly significant at higher values of the cruising speed. 

For example, at a cruising speed of 19 km/h, increasing the traction value from 
1.6x1 O’ 4 to 2.4x1 O' 4 kW h/kg km doubles the mobility power (from 70 to 140 
kW e ). At 30 km/h the same increase in surface traction yields increases the 
mobility power from 153 to 526 kW e (more than a 360% increase). At the lower 
value of 8.0x1 O' 5 kW h/kg km, the mobility power required at 30 km/h is 49 kW e . 
For the 100 kW e base design, a maximum cruising speed of 30 + , 19, and 12.5 
km/h can be obtained for surface traction factors of 8.0x1 O' 5 , 1.6x1 O' 4 , and 
2.4x1 0* 4 kWh/kgkm, respectively. 

3.3 Effect of Reactor System Specific Mass 

The power system specific mass (SM r ) of 100 kg/kW 6i which is used in the 
base case, yields a power system total mass of 10,000 kg. This includes 3000 
kg for the power system based on SP-100 technology. The remaining 7000 kg 
is then used for the radiation shield mass. This shield mass is required to 
insure the crew's safety with a dose rate from the reactor of 30 rem/y during 
operation (justification of this value is given in section 4.1) For the purpose of 
assessing the sensitivity of the total rover mass to the power system specific 
mass over the expected range of possible operation, the specific mass of the 
reactor power systems was varied from 50 to 150 kg/kW e . A more accurate 
estimate of the reactor power system specific mass will be presented in 
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FIGURE 18: Effect of the Soil Traction Parameter on the Mobility Power 





Chapters 4 though 6, when various power plant subsystems are better defined. 
This next phase of the project will focus on the design and integration of the 
radiation shield and the power system with the rover. 

In Figure 19, the total power requirements as a function of cruising speed 
is shown for three different specific mass values (50, 100, and 150 kg/kW e ). For 
all three cases, the minimum power is 30 kW e (user net power only) at 0 km/h. 
At low cruising speeds (< 10 km/h) the power system's specific mass has 
virtually no effect on the reactor power. However, once the speed begins to 
increase beyond 10 km/h, the specific mass begins to show a significant effect 
on the total reactor power requirement. For example, at 30 km/h the total power 
requirements double from 183 kW e at 100 kg/kW e to 378 kW e at 150 kg/kW e . 
Figure 19 also shows that for the base case of 100 kW e , the maximum cruising 
speed attainable for the rover is 23, 19, and 16.5 km/h, for nuclear reactor 
system specific masses of 50, 100, and 150 kg/kW e , respectively. 

Figure 20 shows that the specific mass of the reactor power system; 
however, greatly impacts the total rover mass. At 100 kW e , increasing the 
system’s specific mass from 100 to 150 kg/kW e only decreases the maximum 
cruising speed by 2.5 km/h, but increases the total rover mass from 26 to 32 
metric tons. This mass increase is significant and represents an additional 
launch cost of approximately six billion dollars (at the cost of $1M/kg). 

3.4 Effect of the Utility Cars Mass 

The effect of varying the total mass of the three utility cars on the total rover 
mass and the system power was evaluated for a total RUC mass range from 
15,000 to 18,000 kg. Figure 21 presents a plot of rover mass versus 
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CRUISING SPEED (km/hr) 

FIGURE 19: Effect of Reactor Specific Mass on Total Power Requirement 







•00008 



56 


CRUISING SPEED (knyV) 

FIGURE 20: Effect of Reactor Specific Mass on Total Rover Mass 
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cruising speed for four values of RUC mass. At a speed of 0 km/h the total rover 
mass varies from 18,412 to 21 ,412 kg. This difference in total mass simply 
reflects the increase in RUC mass from 15,000 to 18,000 kg. At 19 km/h the 
same increase in the RUC mass increases the total rover mass from 26,004 to 
30,51 1 kg, a difference of 4,507 kg. The additional 1 ,507 kg reflects the small 
increase in the PPS mass needed to meet the additional mobility power 
requirements. It is therefore concluded that a 300 kg variation in RUC mass is 
not critical in evaluating the PPS power or mass requirements. 

3.5. SUMMARY OF ROVER PARAMETRIC ANALYSIS 

At this point the concept of the manned Mars rover has been developed as 
the layout of a four car train including an auxiliary power system. As was stated 
earlier, the mass of the rover is a key parameter in determining the mobility 
power requirement, which in turn is used for determining the reactor power 
level. Since the total mass of the rover is strongly related to the mission 
scenario, it was necessary to define a reference mission prior to the start of the 
design. This mission, including an emergency return scenario, has been 
developed using the initial design goals stated in section 2 and the mass of the 
individual cars is given in Tables 1 , 4, and 5 for the PCV, EU, and SC, 
respectively. 

In summary, the parametric analysis showed that the total power required 
from a nuclear reactor primary power system and the total rover mass are strong 
functions of the cruising speed, specific mass of the power system, and the 
surface traction and drive efficiency, and to a lesser extent the mass of the rover 
utility cars. The use of an SP-100 reactor, generating 100 kW e , would provide 
for a cruising speed of 16.5 km/h at a reactor system’s specific mass of 150 
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kg/kW e and 19 km/h if the system specific mass could be reduced to 100 
kg/kW q . At this power level the mass of the PPS is 32% of the total rover mass 
for 100 kg/kW e and 53% at a specific mass of 150 kg/kW e . The total mass of 
the rover with a 1 00 kW e PPS and a specific mass of 1 00 kg/kW e is 
approximately 26 tonnes. However, at 150 kg/kW e the rover mass will by 23% 
higher at 32 tonnes. 

With the mass of the PCV, EU, SC determined, the mobility power may 
now be addressed. In the next section, four components of the PPS will be 
studied in order to evaluate the system specific mass. These components 
include the radiation shield, the energy conversion system, the reactor, and the 
waste heat rejection system. The design goal is to integrate the requirements 
for designing the reactor power system with a nominal specific mass of 100 
kg/kW e . At a output power level of 100 kW e the total power system mass will be 
10 tonnes, of which 3 tonnes for the masses of the core, energy conversion, 
heat rejection, and additional subsystems and 7,000 kg in for the biological 
shield. 
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4. RADIATION SHIELD DESIGN 


4.1 Shield Overview 

The radiation shield for nuclear reactor, that is to be operated in the vicinity 
of a human crew, will makeup the largest portion of the overall system mass. 
Many factors must be considered in the design of any radiation shield, some of 
these factors include: 

- size and nature of the nuclear power source, 

- type of radiation produced a, p, y, neutrons, 

- configuration of the spacecraft and payload, 

- length of mission, and 

- total radiation dose limit that will be permitted. 

Although many of the values for these unknowns have already been 
determined by the original mission requirements, several questions remain 
unanswered. For instance, it is already known that the system is to be a reactor; 
therefore, the primary shielding requirements will be for y and neutrons. While 
the configuration of the rover has been developed, much leeway in the actual 
shield configuration and reactor car layout still exists. While a shadow shield is 
preferred for space applications, it is not as appropriate for a Mars rover . In 
space, vacuum conditions eliminate any significant scattering source; however, 
both the soil and atmosphere of Mars will contribute to an increased radiation 
dose rate at the crew's location if precautions are not taken. By increasing the 
shield angle and placing at least a moderate shield around the entire reactor, 
significant dose to the astronauts and prohibitive constraints on crew activities 
will be avoided. 


60 



The first step in designing the radiation shield is to determine the maximum 
allowable dose to the crew. The maximum permissible dose for the astronauts 
is given by the NCRP as 0.5 Sv/y (50 rem/y) from combined manmade and 
natural sources (i.e.cosmic radiation and solar particle events). For a three year 
Mars mission (1 year on the surface) a total received dose for the mission of 1 .5 
Sv (150 rem) can be received. The breakdown of radiation from natural 
sources can be seen in Table 7 [Bloomfield, 1990]. The highest annual 
expected dose rate of 0.2 Sv/y, while on the planet surface, is of principle value 
concern for the manned rover. Using the NCRP's recommendation that the 
maximum annual dose limit be no more than 0.5 Sv/y for astronauts 
undertaking exploratory missions, the dose obtained from all nuclear sources 
onboard the rover should not exceed 0.3 Sv/y (30 rem/y). 

It is important to note that 0.2 Sv/y is the highest expected rate on the 
surface, while 0.1 Sv/y is the average expected rate on the surface. The dose to 
the crew from natural radiation sources while on the surface of the planet is 
significantly less than that received during transit to Mars. This is due to the 
shielding of cosmic and solar particle radiation by the C0 2 atmosphere (~10 
g/cm 2 ). From this information a maximum dose rate was set at 0.3 Sv/y (0.032 
mSv/h). While this dose rate coupled with the maximum surface dose of 0.2 
Sv/y from natural sources gives a total rate of 0.5 Sv/y, this is again the 
maximum value expected. The effect of the dose rate on the required thickness 
of various shield materials and shield mass will be the subject of discussion 
later in the report. 

4.2 Space Shielding History 

Shielding designs have been extensively studied in the late 60’s and early 
70's during both the SNAP and rover Programs. With the appearance of the 
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SP-100 space nuclear power program in the early 1980s, again lightweight 
shields for space systems have become of interest. 

4.2.1 SNAP-1 0A 

To date, the only flight experience with a space reactor shield has been on 
the SNAP-1 OA. This system was launched on April 3, 1965. The reactor used a 
98 kg shadow shield consisting of a cold-pressed Lithium Hydride material 
which was reinforced with a stainless steel honeycomb contained in a SS-316 
casing [Angelo and Buden, 1985]. The design of the shield was to limit the 
radiation dose at the spacecraft interface to 5x10 12 nvt fast neutrons and 4x1 0 7 
rads of gamma radiation (the actual measured dose was significantly higher, 
apparently due to a higher than expected neutron flux caused by scattering from 
the control drums). The operating temperature of the shield was bounded 
between 712 and 579 K. However, the mission was experimental in nature and 
the shield was not designed to protect any human crew. 

4.2.2 ORNL Heat Pipe Reactor Shield 

Since the time of SNAP-1 OA many studies into the shielding for space 
reactor systems have been completed. One such study was performed at Oak 
Ridge National Laboratory (ORNL) for their stainless steel, Potassium Heat Pipe 
Reactor concept. This shield design, shown in Figure 22, uses a combination of 
tungsten (for gamma shielding) and lithium hydride (for neutron shielding). In 
addition, it was shown that to achieve the lowest dose rates the tungsten should 
be placed in two layers separated by a layer of lithium hydride. Due to 
tungsten's large (n.y) cross section, it is important to limit the neutron flux within 
the tungsten in order to minimize production of secondary gammas. According 
to the ORNL report, changing the dose rate requires a change only in the 
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Figure 22: ORNL Heat Pipe Reactor 
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thickness of the first layers of tungsten and lithium hydride. The effect of various 
shield material thickness on the dose rate is shown in Figure 23 
[Barratino.1985]. The dose constraints for the ORNL shield were .03 mSv/h 
within the shields shadow (45°) and 3.0 mSv/h outside the shadow at a 
distance of 32.81m [Engle et al., 1971). 

The Oak Ridge shield study consisted of optimization analysis of both two 
and three cycle shields (a cycle refers to the combination of a tungsten and 
lithium hydride layer). The optimization of the shield was for a 45° zone of 
protection for the two cycle shield and a 90° cone angle for the three cycle 
shield. In order to perform the optimization calculations the ASOP code (ANISN 
Shield Optimization Program) was used. 

For the two cycle shield, the reactor power was taken at 450 kW t with dose 
constraints at the crew quarters of 0.0075, 0.03, and 0.12 mSv/h. The crew 
distance for the design was 61 m (200 ft) from the core centerline. The 
thickness of tungsten and lithium hydride versus the dose rate at the crew cabin 
is shown in Figure 24 [Engle et al., 1971]. There are two interesting points 
about this curve. First, in order to decrease the dose rate it is necessary to 
increase the thickness of the first layers of tungsten and LiH, while keeping the 
second layer of tungsten -1 1 cm and of LiH -42 cm. Next the total thickness of 
LiH required to adequately protect the crew is -70.5 cm, which turns out to be 
about the same thickness as needed to protect the electronic equipment for the 
SP-100 reactor [Barratino,1985]. 
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Figure 23: Shield Material Thickness vs Dose Rate 
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Figure 24: Thickness of W and UH Layers in a Two-Cycle Shield as a 
Function of Dose Rate at a Distance of 200 ft (450 kW t ) 
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4.3 Shield Material Selection 

One of the main objectives in any space system design is the minimization 
of mass. Of course the farther and faster one needs to travel the more costly it 
becomes to launch large heavy objects. Since man-rated shields are by 
definition heavy, every effort should be made to minimize its mass. 

One very attractive alternative to transporting large quantities of material 
from earth, is the use of indigenous planetary soil materials for the shield. For 
Mars the estimated soil composition, as relayed by the Viking landers, is shown 
in Table 7 [Bloomfield]. This table shows 50% oxygen content in the soil. The 
oxygen is primarily found in the form of oxides. Because these estimates are 
based solely on the Mars Viking landing sites, considerable variation may be 
found at different surface locations. Also, the Viking probes only scratched the 
surface in their probing of the Martian soil; therefore, fluctuation may be found 
by excavating down several meters. Although the use of indigenous planetary 

Table 7: Martian Soil Composition 


Element 

Percentage 

0 

50 

Si 

21 

Fe 

13 

Mg 

5 

Al 

3 

S 

3 

Cl 

1 

Ti 

0.5 

K 

0.2 

V 

70ppm 

Sr 

60ppm 

Rb 

30ppm 

Zr 

30ppm 
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surface materials to provide radiation attenuation has been shown to be a 
viable shielding alternative for both Martian and Lunar base power systems, 
showing significant reduction in the required mass to be transported from earth 
[Bloomfield], it is not as attractive for a manned rover shield. 

The minimum thickness of soil material needed to attenuate the radiation 
to acceptable levels would be ~7 m. This, in turn, would have an equivalent 
mass in the 10s of tonnes. While building such a shield could be accomplished 
with modest material handling equipment (which should be available from other 
planetary surface base construction) [Bloomfield, 1990], such a massive shield 
would severely inhibit the movement of a manned surface rover. In addition to 
the increased mass of the operating rover, the transmutation of the soil by 
neutrons produced in the reactor would cause a degree of activation of the 
shield. Table 8 shows the activation of the soil (in curies) for given time after 
reactor shutdown, after 10 years of reactor operation. It should be noted that 
significant radionuclide buildup is present. Although the rover power system 
would not operate continuously for 10 years, the activation of the soil would still 
cause an additional radiation dose to be received by the crew. For these 
reasons it was decided that a rover designed with a shield constructed on earth 
would best suit the mission. 


Table 8: Nuclide Buildup 


Radionuclide 

1 00 days 

400 days 

Radionuclide 

1 00 days 

400 days 


curies 

curies 


curies 

curies 

P-32 

587 

592 

Mn-56 

37.9 

37.9 

AI-28 

361 

361 

Si-31 

87.1 

87.1 

Mg-27 

118 

118 

Ar-37 

43.3 

50.3 

Fe-55 

7.9 

28.3 

Na-24 

27.4 

27.4 

Mn-54 

21.1 

62 

All others 

34.2 

92 
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4.4 Calculation Method 

In order to model the core/shield configuration, the TWODANT neutral 
particle transport code developed at Los Alamos National Laboratories (LANL), 
was used [O'Dell et. al, 1989]. The first step in the shield design was to 
adequately model the reactor core so that an accurate source term could be 
developed for the shielding calculations. Since the core technology to be used 
is SP-100 derivative the current SP-100 design was felt to be the best possible 
choice for the reference model. In order to develop an adequate representation 
of the SP-100 core in a TWODANT input file a model of the core developed by 
John McGhee for helium production analysis was modified and used [McGhee, 
1989]. The material zones are shown in Figure 25. 

4.4.1 Neutronics Model 

Since the model input file was originally developed to be run on a Cray 
computer at Los Alamos, the regional zones had to be simplified so that it could 
be run using the capabilities of a VAX-8750. The model is a two-dimensional r- 
z cylindrical representation of the SP-100 reactor core. Further modification 
was made to the core by replacing the beryllium reflector shutter panels by 
rotating BeO/B 4 C control drums. This was done because, unlike the SP-100 
reference mission which employs a shadow shield in a space environment, the 
rover mission will require a minimum of an oriented 4-rc shield to insure 
adequate protection for the astronauts. Such a shield would severely limit the 
usefulness of a shutter control system. 

4.4.2 Model Verification 

In order to insure lhat the changes to the core model would have a 
minimum effect on the energy distributions, the core model was compared to the 
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more detailed original model. The comparison of the two models is shown in 
Figure 26 (a&b). From the neutron energy spectrum, and Figures 27 (a &b) and 
28(a&b) ( the radial and axial flux distributions) it can be seen that the simplified 
model corresponds well with the results obtained by McGhee, although the 
model used a 42 group cross section library instead of the 80 neutron group 
library originally used. The 42-group library is a coupled neutron (30 energy 
groups) and gamma (12 energy groups), which allows the the gamma spectrum 
in addition to the neutron distributions to be plotted. Figure 29 plots the 
normalized gamma energy distribution at both the core centerline and the 
reflectors outer edge. As is expected the energy distribution, is shifted toward 
the lower end of the spectrum at the edge of the reflector. Figures 30 and 31 
show the normalized radial and axial gamma flux distributions, respectively. 
Although there is no such distribution calculated for the more detailed model ( it 
uses only a neutron cross section library), these curves can still be used to 
evaluate the reasonability of the model. 
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Figure 26b: Neutron Energy Distribution for Original SP-100 Core 



Figure 26a: Neutron Energy Distribution for Modified SP-100 Core 
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Figure 27b: Radial Neutron Flux Distribution for Original SP-100 Core 



Figure 27a: Radial Neutron Flux Distribution for Modified SP-100 Core 
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Figure 28b: Axial Neutron Flux Distribution for Original SP-100 Core 



Figure 28a: Axial Neutron Flux Distribution for Modified SP-100 Core 
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Figure 29: Normalized Gamma Energy Distribution for Modified SP-100 Core 
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Figure 31 : Axial Gamma Flux Distribution for Modified SP-100 Core 
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4.4.3 Dose Rate Calculation Method 

Before going into the analysis of the shield, it is important to discuss the 
manner in which the dose rate is calculated. The output from TWODANT is in 
the form of fluxes at the various energy levels related to the energy groups in 
the cross section library. It is, therefore, necessary to convert these fluxes into 
equivalent dose rates at the location of interest. The calculation of the 
contribution to the dose rate from gammas and neutrons is carried out in two 
different manners. 

4.4.3. 1 Gamma Dose Rate 

Output from TWODANT model of the SP-100 core is in the form of neutron 
and gamma fluxes at specified locations for each of the various energy levels 
related to the 42 energy groups in the cross section library. These fluxes are 
then converted into equivalent dose rates at the location of interest. For the 
gamma dose received the following equation [Lamarsh, 1983] was used: 

D= 5.76x1 0 1 4 4> Y E Y (p a /p)tisQ mSv/h, (11) 

where <J> y is the gamma flux over a particular energy range, E y is the gamma 
energy, (p a /p) t j S is the mass absorption coefficient for tissue, and Q is the quality 
factor for gamma radiation (=1 ). The quantities for (p. a /p) t j S are given in 
Appendix C. 

4.4. 3.2 Neutron Dose Rate 

Dose computations for neutrons are more difficult than for gamma rays due 
to the complex way in which neutrons interact with matter. Because of this 
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Figure 32: Neutron Flux for an Equivalent 0.01 mSv/h Dose Rate 





complex interaction, the energy deposition is most conveniently computed 
numerically using Monte Carlo Method [Lamarsh, 1983]. From this type of 
calculation the curve in Figure 32 can be developed. This curve gives the 
neutron flux ($') as a function of the neutron energy required to deliver a dose of 
0.01 mSv/h. Figure 32 can then be used to compute the dose rate from a given 
neutron flux. To calculate the dose, the curve in Figure 32 was divided into 
three regions. For neutrons with energies in either region 1 and 3 the flux 
needed to produce 0.01 mSv/h was relatively constant at 4>’ * 270 and <>' = 6.8 
n/cm 2 sec, respectively. Region 2 is approximated by the following polynomial: 

<|>'= exp[l .0884 - 0.28082(3 - 0.2256(3 2 - 0.0307 1(3 3 - 0.01 31 p 4 ], (12a) 

and the neutron dose is given by, 

D n (mSv/h) » (0n/<Dn)(1 00). (12b) 

Using equations 11 and 12(a and b), along with the constant neutron fluxes for 
regions 1 and 2, a fortran program was developed to calculate the gamma and 
neutron components of the dose rate from the given fluxes generated by the 
TWODANT code. 

4.5 Shield Design Approach 

Now that the core model has been modified and a method for converting 
the fluxes into an equivalent dose rate has been developed it, is time to 
determine the shield geometry and composition. For the initial evaluation 
lithium hydride was used for neutron attenuation, while tungsten was used to 
reduce the gamma component of the dose rate. 
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Figure 33: One-Dimensional Shield Geometry 



4.5.1 One-Dimensional Layout 

The approach for the shielding design consisted of several distinct steps. 
First, due to the size and complexity of the two-dimension mesh, a majority of 
the cases were run using a one-dimensional model. The results were then 
adjusted to fit the 2-D cases. The one-dimensional geometry is shown in Figure 
33. For all the cases a 3 cm thickness of Al and 1 0 cm thickness of water were 
placed between the reactor and the astronauts. The aluminum was used to 
simulate the walls of the rover and the equipment stored between the reactor 
and the astronauts (this includes the energy conversion equipment which will 
be discussed latter in the report) while the water represents the storage of 
consumable and waste water stored in the rover. The location or the Al was 
12.5m from the core centerline and the water was placed directly in front of the 
Al. The 12.5m location was chosen because it is the midpoint between the 
reactor and the average distance of the crew. The region from the edge of the 
shield to the Al is filled with C0 2 to simulate the martian atmosphere. While 
from the edge of the water to the crew (25m) is air which simulates the 
atmosphere within the pressurized cars (see Figure 3). 

4.5.2 Lithium Hydride Thickness 

Using the one-dimensional model, the first step that was performed was to 
determine the effect of altering the thickness of lithium hydride. This was 
accomplished by first setting the tungsten thickness to zero. The shield radius 
(R shield) was then varied from 50-150 cm. The results from these runs are 
displayed in Figure 34. Both the dose rate from neutrons and gammas can be 
seen in Figure 34. While both components of the dose rate decrease 
proportionally to the log of the LiH thickness, the neutron portion decreases at a 
significantly greater rate. A thickness of 125 cm was chosen to be held constant 
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Rgure 34: Dose Rate verses Lithium Hydride Thickness 




for judging the effect of varying thickness of tungsten on the dose rate. 1 25 cm 
was chosen as the LiH thickness because it relates to a dose rate at the 
astronauts of 1/10 the total allowed rate of 0.032 mSv/h. Since the density of 
tungsten is 19.3g/cm 3 as compared to the 0.73g/cm 3 for lithium hydride, it is 
preferred to use LiH over tungsten whenever possible. An equation will be 
developed, in a later section, which can be used to determine the required 
thickness of LiH once the amount of tungsten has been determined. 

4.5.3 Tungsten Thickness 

By holding the LiH thickness constant the effect on both the neutron and 
gamma dose rates of increasing the thickness of the tungsten layer are studied. 
The results from varying the W layer from 0 to 40 cm are shown in Figure 35. 
The dose rate due to the neutrons shows a dramatic decrease, dropping six 
orders of magnitude. The decrease is exponential in nature which is what is 
expected since W-182 (26.3% W nat ) is a strong neutron absorber. The cross 
sectional data for tungsten is located in Appendix E. 

While the gamma dose also decreases exponentially, this decrease is 
much less pronounced than for the neutrons (three orders of magnitude 
compared with the six for neutrons). This reduced effectiveness for y 
attenuation is not due to W inability to stop the gammas produced by the core, 
but rather the large neutron cross section of the tungsten and secondary 
gammas produced in the shield by the (n,y) reaction. 

4.5.4 Effect of Separation Distance between Tungsten Layer and Core 

In order to reduce the level of secondary gammas produced within the 
tungsten it is necessary to reduce the number of neutrons which reach the 
tungsten. The first step in doing this was to locate the tungsten further from the 
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Figure 35: Dose Rate verses Tungsten Thickness 
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Figure 36: Dose Rate verses Tungsten Distance from the Core 




reactor and inserting a layer of LiH within the gap. Doing this has the 
advantage of decreasing the neutrons by neutron attenuation within the LiH. By 
increasing the distance from the core edge, as expected, the gamma 
component of the dose rate was dramatically reduced. Placing 45 of the 125 
cm of lithium hydride between the reactor and the tungsten lowered the gamma 
dose rate by over four orders of magnitude independent of the tungsten 
thickness. This result is shown in Figure 36 for a tungsten thickness of 10 cm, 

20 cm, and 40 cm. While the use of 40 cm of tungsten placed directly adjacent 
to this reactor resulted in a gamma dose of ~ 0.045 Sv/h (this would surpass the 
0.3 Sv/y limit in < 8 h and would prove fatal to the crew in — 9 days), a distance 
of 45 cm was enough to bring the dose down below the limit of 0.03 mSv/h. In 
fact, at this separation distance, a thickness of just over 20 cm of tungsten would 
be sufficient for a gamma shield. The problem here is that while placing the 
tungsten further from the core reduces the thickness needed, the mass of the 
tungsten for such a shield would be over 30 tonnes. 

4.5.5 Multilayers of Tungsten 

The next step in reducing the shield mass was to divide the tungsten 
gamma shield into 2 layers, similar to the ORNL shield discussed earlier. By 
doing this, the mass of the tungsten may be reduced further. In this case the first 
layer of tungsten was placed adjacent to the BeO reflector outside the reactor 
(22 cm from the core centerline). The separation distance between the first and 
second tungsten layers (£ in Figure 33) was then varied from 0 to 45 cm. In all 
cases the total LiH thickness was 125 cm and the total tungsten thickness 
(xi +X 2 ) was 30 cm. The results of four cases are shown in Figure 37. These 
four cases represent a second tungsten layer thickness (X 2 ) equal to 30 cm, 20 
cm, 10 cm, and 5 cm. The curves of Figure 37 show some interesting 
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Figure 37: Dose Rate verses Tungsten Layer Seperation Distance 



characteristic behaviors of the tungsten portion of the shield. First at £ = 0 all 
four curves converge to the same dose rate. This is expected since the point 
simply represents a single layer 30 cm thick layer of tungsten placed against the 
BeO reflector. 

As the separation distance begins to increase the three split tungsten 
shields show much smaller dose rates than single 30 cm layer of tungsten. The 
smaller dose rate is due to the fact that the first layer of tungsten, while 
producing many secondary gammas, also removes (through absorption) many 
of the neutrons produced in the core. In addition many more neutrons are 
removed by the LiH section separating the two W layers. The second layer of 
tungsten then absorbs many of the secondary gammas while producing few 
additional secondary gammas due to the lower neutron flux. 

As the separation distance (£) between the W layers is increased beyond 
10 cm the curve for X2=5 cm levels out at a constant 1 .5 mSv/h while beyond 30 
cm the X2=10 cm curve levels out at 0.003 mSv/h. This is because as £ is 
increased a point is reached where the neutron flux is at a point where 
secondary gammas produced in the second layer becomes negligible. When 
the thickness of the second layer becomes larger, the separation distance at 
which the neutron interaction becomes negligible also increases. This is not 
only because the neutrons have a higher probability to interact within the larger 
thickness, but more importantly because the first layer is smaller; therefore, the 
flux is greater at the front of the second layer of tungsten. As seen in Figure 37, 
the combination of a 20 cm first layer and a 10 cm second layer gives the lowest 
dose rate up to % =38 cm where it crosses the X 2=20 cm line. This is preferred 
because not only is the dose rate lower, but the mass is also lower since a 
larger fraction of the tungsten is at a smaller radius. 
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Figure 




4.5.6 Dose Rate Results for Reduced Tungsten Thickness 

Since the dose rate for the X 2=10 cm shield is 8.5 times smaller than the 
maximum dose rate, the shield mass can be reduced by eliminating a portion of 
the tungsten. In Figure 39 data, similar to that in Figure 38 is plotted; however, 
in this case, the sum of both tungsten layers is 20 cm. Here the dose rate to the 
astronauts is 0.03 mSv/h due to gammas when x-| =10 cm, X2=10 cm, and £ =35 
cm. Using these dimensions will yield a substantial weight savings over the 
previous 30 cm thick tungsten shield. 

4.5.7 Multilayer Tungsten/ Depleted Uranium Shield 

The final step in order to evaluate a minimum shield mass is to replace the 
second layer of tungsten with depleted uranium (U-238). In this stage, the 
second layer of tungsten in the x-) =20 cm X 2=10 cm shield was replaced by 
depleted uranium. A layer of U-238 is chosen to compare to tungsten because 
of its excellent gamma attenuating abilities [Lamarsh, 1983] and its lower 
density (18.5 g/cm 3 compared to 19.3 g/cm 3 )[Walker, et.al, 1984]. The 
disadvantage with U-238 is that if it is exposed to fast neutrons (E n > 0.6 MeV), 
the atom can undergo fission, causing an increased dose rate to the astronauts 
by both secondary gammas and neutrons. In addition, the fissioning of U-238 
will increase the heat generated in the radiation shield, which may cause 
severe deterioration the shields ability to attenuate the oncoming radiation may 
occur if the shields temperature becomes too high [Barratino, 1985]. For these 
two reasons only the second layer of the shield used U-238. Being placed 
outside a layers of tungsten and two layers of lithium hydride the neutron flux 
(especially fast neutrons) is significantly lower. 

When the second layer of tungsten is replaced by depleted uranium the W- 
W shield showed a slightly lower dose rate. The results of the W-U shield 
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Rgure 39: Comparison of Astronaut Dose Rate W/W and W/dep. U 
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Figure 40: Two-Dimensional Core Shield Configuration 






analysis are compared with those from the previous W-W shield in Figure 39 (a 
and b). While the curve shows that a combined tungsten-uranium shield is 
superior for a x-j =20 cm and X2=10 cm, when x-j is reduced to 10 cm the W-W 
shield yields a significantly lower dose rate, therefore, a W-W shield was 
chosen for lower shield mass. 

4.6 Two Dimensional Verification 

For the two dimensional verification of the shield layer thicknesses, the 
core shield configuration shown in Figure 40 was used. For all the runs the r 
direction is identical to that used in the one-dimensional case. The core model 
previously described was inserted as the source generator, and a secondary 
layer of tungsten was placed above and below the reactor. The first tungsten 
layers above and below the reactor are smaller that x-j . The justification for this 
is twofold. First, as can be seen by the flux curve in Figures 28 and 30 the 
radiation levels from both neutrons and gammas are lower at the ends than at 
the sides. Secondly, a higher dose rate at the edge of the shields may be 
allowed, with only sufficient shielding needed to protect the astronauts from 
radiation scattered from the CO 2 atmosphere or the Martian soil. The results of 
the two-dimensional runs are shown against the 1-D results in Figure 41. It can 
be seen that going to the two dimensional geometry reduces the calculated 
dose to the astronauts by ~15%. This helps to confirm that the 1-D dimensions 
calculated for the shield are more than sufficient to protect the astronauts. 

4.7 Shield Design 

From the results in the previous sections it can bee seen that a cylindrical 
shield consisting of 63 cm of lithium hydride, and two layers of tungsten, the first 
being 1 1 cm thick while the second is 9 cm thick will sufficiently protect a group 


94 



2-D Curves 



0 5 10 15 20 25 30 35 40 45 

Separation Distance (cm) 

Figure 41 : Comparison of 1-d to 2-D Results 
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of astronauts at a distance of 25m. However, while this shield reduces the 
radiation from the reactor to a sufficiently low level for the crew, it is be 
extremely massive and protects a larger area than needed. 

For the train rover it is anticipated that a majority of the astronaut's time will 
be spent within the rover. Therefore, it is of primary concern to protect this area. 
If the shield is designed so that some of its mass on the sides not directly facing 
the astronauts is trimmed away, then a large mass saving may be incurred. 
Because of the vacuum of space removing large portions of a shield (leaving 
only a shadow of protection) mearly limits the area which is usable by 
equipment and crew. However, on Mars scattering off of both the soil and 
atmosphere can increase the dose rate. For this reason, along with the need for 
limited access around the reactor car, the shielding on the side of the reactor 
cannot be completely eliminated. 

In addition to the thickness of the side shielding is the forward angle which 
will get the maximum protection from the shield. Figure 42 shows the layout of 
the rover along with the 25m radius from the reactor and several forward 
angles. The minimum shield angle is determined by the smallest possible 
angle needed to completely cover the PCV and EU. Of course, if any significant 
amount of work is to be completed outside the rover, then a larger forward angle 
will be required. 

Besides protecting the astronauts during any work outside the rover, the 
vehicle must be protected while turning, Figure 43 shows the turning protection. 
With the minimum angle (8.6%) primary control vehicle leaves the zone of 
maximum protection (2MP) with only a few degree turn. While a 90° forward 
angle on the shield would give protection for even the most severe turns, it is felt 
that such a small fraction of the transit time would be dedicated to high degree 
turn that this would be too conservative. In addition, by taking proper 
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Figure 42: Area of Protection For Various Shield Forward Angles 




Figure 43: Turn Protection for Various Shield Forward Angles. 




precaution, a majority of the severe turns may be avoided. A forward shield 
angle of 50° was selected. This will allow the rover to have a turning radius of 
40 m while remaining in the ZMP, although sharper turns should be avoided if 
possible, they can be made as long as their time of execution is kept small. 

With the information in the previous curve, the shield configuration and 
dimensions can now be determined. From Figures 34 and 35 an estimate of the 
required lithium hydride thickness may be determined given the total tungsten 
thickness used. The equation is as follows: 

D.R.(LiH) = C exp(T w ) expOm) (13a) 


or solving for the LiH thickness and rearranging, 

T [ D.R.(LiH) 

L,H = [Cexp(T w ). • 


(13b) 


Of course in order to use these equations for the lithium hydride thickness, it is 
first required that the tungsten thickness is known. Also, equation 1 3b is no 
longer valid if the separation distance t, is less than the LiH thickness. 

To determine the required W thickness and the separation distance, £ , 
curves similar to Figures 35 and 36 were used. The typical shield layout is 
shown in Figure 44 and 45. For the base case of a 500kW t reactor we have 
xi *1 1 cm, X2=1 0 cm, £ =40 cm, and T LiH =65 cm for the forward direction. For 
the sides dimensions are slightly smaller due to the higher allowed flux at the 
surface of the shield. Here the dimensions are: 

x 1 = 7 cm 

x 2 = 9 cm 

% m 30 cm 

T LiH -41cm 
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Figure 44: Reference Shield Layout for Manned Mars Rover (top view) 
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Figure 45: Reference Shield Layout for Manned Mars Rover (side view) 




The total shield mass is 16,567 kg. In order to develop an approximation for the 
shield mass, an analysis was performed at four additional power levels (100, 
200, 300, and 1000 kW t ). The mass versus power level is shown in Figure 46. 
By fitting an equation to the curve the following formula was generated: 

SMghield * 1494 (P th )*°' 6203 (14) 

where SM S hield is the shield specific mass in kg/kW t and P th is the reactor 
thermal power in kW t . While this equation seems rather simplistic, it should be 
remembered that shield are quite reactor specific and much more than a rough 
estimate can give a false sense of finality to a conceptual design. By comparing 
masses calculated by this equation to designed man-rated shield masses 
[Bloomfield et. al,1989; Engle et. al, 1971], the results appear to be reasonable 
(in a general sense) to within 1000 kg. This equation will be used later once the 
reactor thermal power is determined from the energy conversion system 
efficiency. 

A source of conservatism, that will be seen in the results of the total PPS 
mass, is the assumption that the reactor will be operating at its nominal thermal 
power all the time and assuming constant core dimensions, regardless of the 
reactor thermal power level. In reality, it is expected that the reactor power at 
night will be reduced by more than 50% when mobility power is not needed (life 
support power and user power are only 30 kWe). When a thermal power of 
0.75 of nominal is used in equation (14) to account for a 50% reduction in 
reactor power during the non-driving periods, the shield mass can be reduced 
by 10 percent without exceeding the maximum permissible dose rate of 0.032 
mSv/h to the crew on board the rover vehicle. 
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5.1 Reactor Specific Mass 

it will be shown that the selection of the conversion system will strongly 
influence the mass the other three main system components; namely, the 
reactor, radiation shield, and waste heat radiator. An equation for the effect of 
the thermal power on the mass of the shield has already been developed; 
however, it is also necessary to have the relationship for the reactor mass 
versus the thermal power of the reactor. Increased power requires a larger 
reactor mass and size. The following equation is developed from a curve fit of 
past and present space reactor designs (see Figure 47) and is used to calculate 
the reactor specific mass as: 

SMreactor = 66-54 (P, h )-°- 7189 (15) 

Although equation 15 is not an exact relationship, it does give a 
reasonable estimate of the reactor specific mass. The maximum deviation from 
the mass of the reactor system used was 15%, which yields an error in the 
calculated mass of 50- 150 kg. This deviation is only a small fraction of the 
overall reactor system mass; therefore, this difference will yield less than a 2% 
error in the calculation of the total system mass. 

5.2 Radiator Specific Mass 

For three of the four energy conversion systems considered in the next 
section the waste heat is rejected by means a thermal radiator. Both the surface 
area and the mass are very sensitive to the rejection temperature, with higher 
rejection temperature corresponding to lower radiator mass. It will be seen; 


104 



REACTOR THERMAL POWER (kWll 


10000 


1000 


100 


10 



6 




however, in the next section that this high radiator temperature will decrease the 
conversion efficiency which increase the shield mass [Angelo and Buden, 

1985]. The manned rover offers both advantages and limitations for a radiator 

heat rejection system. 

In order to maintain simplicity and allow the reactor to operate continuously 
the area of the radiator panels is limited to the surface area of the reactor car, 
<100 m 2 . While this will restrict the lower temperature end of the energy 
conversion system, it will allow the reactor car to remain a single self contained 
unit, without large awkward fins. The radiator material selector for use in the 
rover should exhibit the following qualities: 

1) high strength to mass ratio, 

2) high thermal conductivity, 

3) compatibility with working fluid or coolant, 

4) compatibility with Mars atmosphere, 

5) high emissivity, and 

6) good structural integrity. 

The material considered for the radiator will eventually be dependent on 
the power conversion system, since the correct choice can only be made once 
the operating temperature is known [Angelo and Buden, 1985]. Although the 
material selection cannot be made without a better understanding of the 
operating conditions currently space radiator designs range in a mass per area 
from 1 to 10 kg/m 2 depending upon the amount of armor included. On the 
surface of Mars micrometers (a major concern in space) are not a concern, 

p 

therefore, a large amount of armor is not necessary. A conservative 5 kg/m 
was selected to be used as the radiator mass. 
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The radiator area can be calculated as a function of the radiator 
temperature and the ambient temperature raised to the 4tb power [Bennett and 
Myers, 1982]. For the radiator calculations the following equation is used: 


A 


o 



^ ( e Trad * oT^n) 


(16) 


where A is the radiator area (m 2 ), P re j is the thermal power rejected (kW t ), o is 
the Stefan-Boltzmann Constant, F is the view factor, e is the radiator surface 
emissivity, a is the radiator absorptivity, and T ra( jand T am are the radiator and 
ambient temperatures, respectively. Using values of radiator materials currently 
available [Angelo and Buden, 1985] and information on Mars [Kaplan, 1988], 
the radiator variable are selected as follows: e = 0.85, a = 0.2, and T am = 250 K. 
The view factor F is taken as 0.8 to allow for reflection off the planet’s surface. 
The Specific Area of the radiator versus the radiator temp is shown in Figure 48, 
where the importance of maintaining as high a temperature as possible can be 
seen. 

The specific mass equations for the radiator, reactor, and shield may now 
be used to evaluate the mass benefits obtained from various energy conversion 
systems. This comparison will be carried out in the following chapter. 
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Figure 48: Radiator Area and Specific Mass verses Surface Temperature 



6. ENERGY CONVERSION SYSTEM 


As we have seen in the last sections, the mass of the radiation shield, 
reactor, and waste heat radiator are all dependent on the thermal power of the 
reactor. The reactor thermal power is controlled by two main factors: 1 ) the 
electrical power required by the rover, and 2) the efficiency of the energy 
conversion system. In order to minimize the shield mass and overall system 
mass, an efficient energy conversion scheme is required. In this section four 
energy conversion systems are assessed for use with an SP-100 type nuclear 
reactor system in a manned mars rover. The conversion systems under 
investigation are: 1) SiGe/GaP Thermoelectrics (TE), 2) C0 2 Open Brayton 
Cycle (OBC), 3) He/Xe Closed Brayton Cycle (CBC), 4) Free Piston Stirling 
Engine (FPSE). Each of these four systems display certain advantages and 
disadvantages which will now be discussed. 

6.1 Thermoelectric Conversion 

Thermoelectric (TE) conversion is a way to directly convert thermal energy 
into electrical energy based upon the Seebeck effect. While this effect was 
originally noticed through the junction of two dissimilar metals, the appearance 
of specially tailored semiconductors has increased the conversion efficiency. 
Modern TE devices use the junction of doped p-type and n-type semiconductor 
[Angrist, 1976]. The electrical output of the TE device depends on the 
temperature drop across the TE, the properties of the semiconductor (electrical 
resistance, thermal conductivity, and Seebeck coefficient), and the specific 
design of the converter element. Figure 49 the layout of an idealized converter 
[Angelo and Buden, 1985]. 
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Figure 49: Idealized Thermoelectric Converter. 
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TE converters have been flown in space by the U.S. since 1961 with the 
launch of the of the Transit 4a and 4b satilites. Since then TEs have amassed 
an impressive flight history, having been flown on all of NASA's most exciting 
missions to the outer edges of the Solar System. TEs have demonstrated both 
a high degree of reliability and predictability throughout their lifetime [Angelo 
and Buden, 1985]. 

The use of TE conversion in the area of space nuclear power has, 
however, been primarily limited to RTG systems. The only flight experience with 
TE conversion on a space reactor was on the SNAP-1 OA system in 1965. The 
TE material used was a Silicon/Germanium (SiGe) mix. While the TEs 
performed as expected, the reactor was accidentally shutdown after only 42 
days of operation. The SP-100 Space Nuclear Power System, currently under 
design by General Electric, will also use TE conversion. While the standard TE 
element used in the past has been the unicouple design (see in Figure 50), the 
SP-100 TEs are designed in the form of multicouple cells. The layout for these 
cells are shown in Figure 51. Unlike the Unicouple which contains a single TE 
element per couple, the Multicouple contains 32 elements per cell. The total 
SP-100 power conversion system design calls for a total of 4320 multicouples 
to be used. 

Because power conversion system of the SP-100 has already been 
designed to deliver 100 kW e , the same design parameters are used for the 
evaluation of TE conversion in the rover power system. These Parameters 
include: T rad = 800K, number of cell = 4320, conversion efficiency = 5.4%, 
reactor thermal power = 2331 kW th , rejected power = 2131 kW th , and converter 
mass * 370 kg [General Electric, 1988]. From these values the component mass 
for thermoelectric conversion on a manned rover may be calculated from the 
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Figure 51 : Multicouple Element For SP-100 



equations developed earlier. The results of these calculations are given in 
Table 9. 


Table 9. Mass Breakdown of TE Conversion System for Manned Rover. 

Reactor Shield Converter Radiator 

563 kg 27600 kg 1300 kg 430 kg 

While Thermoelectric systems have shown a high degree of reliability in 
past use from the component mass breakdown it is seen that the mass of a 
rover powered by a nuclear reactor with a thermoelectric conversion system 
comes nowhere near the design goal of 100 kg/kW e . For this reason a 
Thermoelectric system is ruled out as a choice for the rover energy conversion 
system. 

6.2 C0 2 Open Brayton Cycle 

The use of a C0 2 Open Brayton Cycle (OBC) in a Mars rover eliminates 
the Thermoelectrics problem of large shield mass due to low conversion 
efficiency. By using the atmospheric C0 2 as the working fluid the system has 
the potential to operate at extremely high efficiencies (>50%), due to the low 
ambient temperature (t ave ~230 K). 

The OBC operates according to the P-v and T-s diagrams shown in Figures 
52 a and b. Figure 53 shows a system layout for OBC with a regenerator. The 
cycle consists of intake of atmospheric C0 2 , isentropic compression, isobaric 
heat addition, isentropic expansion, and finally expulsion of the gas back to the 
atmosphere. The efficiency of the a simple OBC can be 
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Figure 52a: P*v Diagram For Open Brayton Cycle 
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Figure 52b: T-s Diagram for Open Brayton Cycle 


116 



117 


Figure 53: System Diagram for Open Brayton Cycle 





calculated as : 


( 17 ) 


n,h ” 1 * <P 2 -pV 1/k) 

where is the thermal efficiency, P-j if the inlet pressure, P 2 is the compressor 
outlet pressure and k is the ratio specific heat at constant volume to specific 
heat at constant pressure. This equation shows that the efficiency is a function 
of the isentropic pressure ratio [Van Wylen and Sonntag, 1985]. The efficiency 
versus pressure ratio is shown in Figure 53. From this curve the the OBC is 
attractive because the high efficiency will yield a lower reactor thermal power 
and lower shield mass. In addition since the OBC is an open system there will 
be no need for a waste heat radiator, again decreasing the mass of the entire 
system. 

For the CO 2 Open Brayton Cycle a simplified axial flow turbine model was 
used for calculating the mass and size of the converter. The axial flow turbine 
model determines the blade length and disk radius for each stage of a turbine 
by calculating (1 ) the energy transfer from the working fluid to the rotating 
blades for each stage of the turbine; (2) the flow area required by the working 
fluid as it exits each stage; and (3) the limiting stage blade speed due to 
material strength considerations. The process is done in an iterative fashion 
due to the interdependence of these calculations. 

The model can accommodate various working fluids (including CO 2 and 
turbine structural materials. A nickel super alloy was used as the turbine 
material along with the CO 2 working fluid. However, CO 2 OBC was rejected, 
because at the low atmospheric pressure on Mars (~10 mbar) the compression 
ratio required to achieve the high efficiencies become large(50 to 100) and the 
the converter becomes the dominant mass component. Table 10 lists the mass 
Breakdown for the OBC. The large converter mass again pushes the system 
specific mass far above the design goal of 100 kg/kW e . 
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Figure 54: Efficiency verses Temperature for Open Brayton Cycle 



Table 10. Mass Breakdown for Open Brayton Cycle. 

Reactor Shield Converter Radiator 

302 kg 11725 kg 10580 kg 0 

6.3 He/Xe Closed Brayton Cycle 

The He/Xe Closed Brayton Cycle is the standard Brayton Cycle discussed 
with use as a dynamic space energy conversion cycle. The P*v and T-s are 
similar to those in Figure 52 a and b for the Open Brayton Cycle; However, 
because it is a closed system, instead of expelling the effluent back to the 
atmosphere, there is an isobaric heat rejection from 4 to 1. The system diagram 
for the CBC is shown in Figure 55. By adding a regenerator to the system the 
overall efficiency may be increased above that of a simple CBC. Figure 56 
shows the T-s diagram for a CBC with regenerator. For a CBC with a 
regenerator the efficiency is given by the following equation: 

■ 1 T^f , * 1, (18) 

Here T 1 /T 3 is the ratio of minimum to maximum temperature. 

Two different size engines were used for the evaluation of the CBC. In 
order to maintain a degree of redundancy each of the two system contained an 
additional conversion unit. The mass breakdown for the systems studied is 
given in Table 11. As in the case of the Thermoelectric system the largest mass 
component is the shield. While the mass of each of the two systems is 
significantly lower than the TE system, both still have a specific mass of almost 
twice the design goal of 100 kg/kW e . 
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Figure 55: System Diagram for Closed Brayton Cycle 






Table 1 1 : Mass Breakdown of CBC Conversion System for Manned Rover 



Reactor 

Shield 

Converter 

Radiator 

3x50 kW e 

325 kg 

12900 kg 

4362 kg 

492 kg 

2x100 kW e 

320 kg 

12507 kg 

5133 kg 

491 kg 


6.4 Free Piston Stirling Engine 

The P*v and T-s diagrams for the Stirling Cycle are shown in Figure 56 a 
and b. Heat is transferred to the the working fluid during a constant volume 
process (2-3) and during a isothermal expansion (3-4), while heat is rejected 
during (4-1) and (1-2). The Free Piston Stirling Engine technology is currently 
being developed under the guidance of NASA Lewis Research Center. The 
FPSE is a reciprocating engine that contains only two moving parts. The 
converter uses gas bearings for low friction and wear, along with a unique 
mechanically balance design for low vibration. Figure 57 shows some of the 
features of the FPSE. The system diagram for a FPSE coupled with an SP-100 
reactor core can be seen in Figure 58. Again multiple converter units would be 
used in the system. The efficiency versus radiator temperature is seen in Figure 
59 for different sized 50 kW e engines. 

For calculating the mass and efficiency for the FPSE an empirical 
correlation of the conversion efficiency and the temperature ratio (TR), specific 
mass of the converter (SM C ), and the electrical Power (P e ) was used. The 
correlation was a third degree polynomial with respect to both TR, SM C , and P e 
determined from a preliminary test conducted at Lewis Research Center 
(LeRC). 

The equation was then used with a SMC of 6-kg/kW e over a range of 
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Figure 56a: P-v Diagram for Ideal Stirling Cycle 
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Figure 57: Free Piston Stirling Engine 


125 


Sdooi nniHin 
ivNomaavoi 



126 


Figure 58: System Diagram for Free Piston Stirling Engine 





temperature ratios from 1 .5 to 3.0. The determined efficiency is in terms of 
fraction of ideal (Carnot) efficiency. From the electric power and specific mass 
of the FPSE, the total conversion system mass may be calculated along with the 
efficiency. Subsequently, the reactor thermal power and the heat rejected from 
the radiator are calculated, and have the mass of the entire primary power 
system may be calculated. The results for both 25 kW e and 50 kW e engines 
are given in Table 12. In each case an additional 50 kW e potential is used for 
redundancy. 


Table 12: Mass Breakdown of FPSE Conversion System for Manned Rover 



Reactor 

Shield 

Converter 

Radiator 

6x25kW e 

314 kg 

12200 kg 

2914 kg 

500 kg 

3x50kW e 

309 kg 

12025 kg 

2693 kg 

476 kg 


Once again the specific mass of the system is significantly higher than the 
design goal. The uses of FPSE technology will; however, allow the lowest 
specific mass (150-160 kg/kW e ) of the four systems studied. 

6.5 Parametric Analysis of FPSE and CBC 

In the previous sections the mass was seen to be lowest for the FPSE and 
CBC systems. Because these masses were calculated using projected values 
for the cycle efficiency, the effect of the efficiency on the overall system mass 
was examined. Here the efficiency of the FPSE was taken to be 29% while that 
of the CBC system was 25%, which represent current design capabilities. For 
both systems the converter mass given in Tables 11 and 12 were used. While 
the mass of each system increased, due to the lower system efficiency, they 
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Figure 59: Cycle Efficiency for FPSE 




both had a lower mass than the C02*0BC or the TE systems. The 3x50 kW e 
FPSE showed the lowest mass of 17.6 tonnes with the 6x25 kW e at 17.8 
tonnes. The CBC systems had a slightly higher calculated mass at 20.1 tonnes 
for the 3x50 kW e system and 20.9 tonnes for the 2x100 kW e system. The rover 
mass breakdown for each of the power system options investigated are shown 
in Table 13. 

6.6 Converter Summary 

None of the four systems considered (TE, OBC, CBC, FPSE) was able to 
meet the design goal, set in Chapter 3, of 100 kg/kW e . At this point the best 
alternative from a mass standpoint appears to be the Free Piston Stirling 
Engine, with a specific mass of 150- 160 kg/kW e . While the Closed Brayton 
Cycle power system has a larger specific mass (185-190 kg/kW e ), it should be 
noted that the technology needed to implement this system is much further 
along in its development. The parametric analysis of the FPSE and CBC 
system showed that even at reduced efficiencies, these systems still required a 
lower mass than the CO 2 -OBC and the TE systems. 

The overall rover mass breakdowns for the conversion systems (excluding 
the OBC) are shown in Figure 60 with the operating parameters for each system 
listed in Table 14. As stated before it is seen that the radiation shield is the 
largest mass component for each system. Although the layout of the reactor car 
will vary from component to component Figure 60 shows the layout for 3x50 
kW e CBC system. From the analysis in Chapters 4,5 and 6, a new range for the 
specific mass of the system was set at 150-190 kg/kW e . The rover mass over 
this range is shown versus the cruising speed in Figure 61. By reducing the 
maximum rover speed, from the base case value of 19 km/h, the mass and 
power requirements of the rover will be kept to within a reasonable level. 
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Table 13: Rover Mass Breakdown for Different Power Conversion System Options 
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Table 14: Operating Parameters and Calculated and Projected Efficiencies 
for Various Energy Conversion Systems. 


Converter Number of 
Units 

Temp. 

Ratio 

Radiator 

Temp. 

K 

Pressure 

Ratio 

Cycle 

Efficiency 

% 

Projected 

Efficiency 

% 

25-kWe FPSE 

6 

2.7 

500 


40 

25 

50-kWe FPSE 

3 

2.7 

500 

- 

42 

27 

50-kWe CBC 

3 

3.16 

395 

1.667 

36.1 

23 

100-kWe CBC 

2 

3.21 

389 

1.662 

37.1 

25 

TE< a > 

4320 

1.5 

800 

- 

>5 

5 

co 2 -obc 

2 

3.5 

- 

2 

43 

28-30 


( a )Number of Thermoelectric Conversion Assemblies (TCA), each contains a 
total of 32 SiGe/GaP Thermoelectric unicouples. 



Figure 61 : Reactor Car Layout for CBC using 3x50 kW e Turboaltinators 
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7. SUMMARY AND CONCLUSIONS 


The study found that the life support power requirement per astronaut 
strongly affects the mass and size of the PV panels for the auxiliary power 
system. Therefore, a careful assessment of life support power requirements is 
strongly recommended. The mass and the size of the PV panels depends on 
parameters such as the emergency cruising speed, the surface density and 
efficiency of the PV panels, and the number of astronauts on board the rover. 
Results show that in order to minimize the size and mass of the PV panels, 
traversing during emergency return should be made at night and day time travel 
should be limited. 

The primary power requirement of 100 kWe was met using an SP-100 
type nuclear reactor coupled to a dynamic conversion system and a man-rated 
radiation shield. The masses of the shield, reactor, and the radiator were 
directly dependent on the reactor power, while this power was dictated by the 
electrical output and the conversion system parameters. The following four 
conversion systems were investigated for use in the PPS: Free Piston Stirling 
Engine (FPSE), He/Xe Closed Brayton Cycle (CBC), C0 2 Open Brayton Cycle 
(OBC), and Thermoelectric conversion. Of the systems studied, the FPSE and 
CBC systems yield the lowest overall PPS mass, with the Stirling system (15-16 
tonnes) giving slightly lower system mass than the Brayton system (18-19.5 
tonnes). Because the shield was such a large portion of the overall system 
mass, its size must be kept as low as possible. Therefore the system with of the 
highest conversion efficiency has the lowest overall mass. This associated 
mass of the 150 to 190 kg/kW e range for specific mass is shown in Figure 61 . 
The overall rover mass for such a system could range from 30 to 34 tonnes. 

Since the efficiencies calculated for both the Stirling and Brayton cycles 
were projected upper limits a parametric analysis was completed to see what 
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effect lower efficiencies would have on the total rover mass. The analysis show 
that while the masses of each system increased by reducing the efficiencies to 
29% for the Stirling and 25% for the Brayton, both systems still had a total mass 
less than the thermoelectric and the open Brayton systems. 

It was seen that major mass component for the PPS was the radiation 
shield. While it was possible to design an effective man-rated radiation shield 
(satisfying the dose rate limit of 0.3 Sv/y to the rover crew) for a manned Mars 
rover using a single layer each of tungsten and lithium hydride, the shield would 
not be mass efficient. A double layered W and LiH shield was more effective 
and significantly lighter. The thickness and positions of the tungsten layers, 
within the shield, were varied to determine the impact of secondary gammas 
produced in the W on the dose rate and total shield mass. 

The combined shield of LiH and W/depleted U was slightly heavier than a 
double layered W-LiH shield. The mass of the optimized man-rated W-LiH 
shield for the Mars rover vehicle powered by a nuclear reactor system increases 
with the thermal power of the reactor raised to the 0.3797 power, and varies 
from 8600 to 20580 kg over the power range of 100 to 1000 kWth, respectively. 
These shield masses were based on a conservative design correlation 
developed based on the results of 1-D neutronics analysis. Results show that 
shield mass estimates could be as much as 15 percent higher than those based 
on 2-D neutronics calculations. A manned Mars rover can easily be designed 
with a SP-100 type reactor power system; however, the astronauts would be 
limited in their freedom by the forward angle of the shield. 
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8. AREAS OF FUTURE WORK 


Since this project was a conceptual design there is still considerable work 
that needs to be done in order to finalize any one for the systems considered. 
Both structural analysis and component integration need to be addressed in 
great detail. In addition much work is still needed in the area of shielding. No 
thermal analysis was performed on the shield and it may be necessary to 
design an active cooling system to keep the maximum shield temperature within 
an acceptable level. 

In this study it was assumed that the power system components could be 
sufficiently isolated from the Martian atmosphere (either by coatings or sealing 
off the reactor car) that corrosion would not affect them over the rover lifetime It 
may be also necessary to investigate the compatibility of any such isolating 
materials and the effect of high temperature upon them. In addition the 
adaptation the Mars rover for lunar applications may show new advantages in 
applicability and economy. 
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APPENDIX A 

MODEL DESCRIPTION AND FLOW CHARTS FOR ROVER MASS PROGRAM 
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MANNED ROVER POWER/MASS MODEL 
Flow Chart 
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SUBROUTINE INPUT 


Denotes Subroutine 
Has not been completed 
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PVFC SUBROUTINE 


Flow Chart 
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ROVOUT SUBROUTINE 
Flow Chart 


Input From 
Main 


Dimension 

Variables 


r Print 
Input Data 


sp - Rover speed. 

SMR - Specific mass of the reactor, 

mnspd - Minimum rover speed for a 
variable speed scenario, 
mxspd * Maximum rover speed for a 
variable speed scenario, 
minspd * Minimum rover Speed for 

multiple single speed scenarios, 
maxspd - Minimum rover Speed for 

multiple single speed scenarios. 


Rover Speed 
Constant? > 


l sp= minspd is \ NO 

sp s maxspdV 
V sp = i+1 / 


Output 

mnspd, 

mxspd 
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ZEROIN SUBROUTINE 
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SPDAVE SUBROUTINE 



tvMpd«0 
totime«0 
minapd»^d(1) 
maxapd « «pd(l) 


I 


r " " 

1*1 


i m N-1 

t 



^TES 


lavespd 

* */espd | 


+«pd(1)*tm(i) | 




1 


sdavspd * aveepd/24.6 
avaepd - avespdlolime 



\ 

If 

i* count 

> 

iK um 

Ives 


.NO 


maxapd * tpd(i) 


mintpd*tpd(i) 







Calculate 
Reactor Power 

ii * — 

■1 


EZ 





*pd() 

- Rover speed over time interval. 

ave&pd 

- Average rover speed. 

savspd 

- Average daily rover speed. 

M) 

- Time intervals. 

totime 

- Total mobility time. 
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Program Reactor Hass 


c 

c 

real rovmas ( 0 : 50 , 5 ) , pa, x, mpcvl, mpcv2, meu, mac, v, 

& rovp(0:50,5) , mpv, mrov, fc, smr(0:50,5), 

t apd(0:50), tm(0:50), mxspd, mnapd 

integer i, j, range, aystyp, modtyp, auxtyp, smtyp, ayat, 

& count, minapd, maxapd 

character*60 xaxia, yaxia, title, conac(ll), intype, 

& rerun 

c 

print*,' Do you what an input prompt (yea or no)?' 
read( 5, 10) intype 

open(unit-40,file-'pvmas' ,status-'new' ) 

10 format(a) 

if ( intype. eq. 'no' ) then 

open (unit-16, file-'rdat' ,atatua-'old' ) 
rewind 16 

read(16,*)pa, amr(l,l), mnapd, mxapd, x, mpcvl, mpcv2, meu, 

& mac, na, modtyp, aystyp, amtyp, auxtyp, ayat 

if (modtyp. eq. 1 ) then 
read( 16 , * ) count 
tmck - 24.6 
do 11 i2 - 1, count 

read( 16 , * ) tm(i2), spd(i2) 
tmck - tmck-tm(i2) 
print*, tm(i2), spd(i2) 

11 continue 

if ( tmck . gt . 0 . 01 ) then 
count - count+1 . 0 

tm( count) - tmck 
spd( count) - 0.0 
endif 
endif 
else 

call rovdat(pa, smr(l,l), mnapd, mxspd, x, mpcvl, mpcv2, meu, 

£ mac, na, modtyp, aystyp, amtyp, auxtyp, ayat, spd, 

£ tm, count) 

endif 
c 

open(unit-35,file-'pwrout' ,status-'new' ) 

j - 1 

15 write ( 35 , 16 ) j 

16 format!/, 35x, 'RUN #' ,i2,///) 


if ( auxtyp. eq. 1 ) then 
f c -0.0 
fcef - 1.0 
mpv - 0.0 

else 

call pvfc(na, mpcvl, mpcv2, x, mpv, meu, mac, fc, fcef) 
endif 

print*, 'fc - ',fcef 

mrov - mpcvl+float(na) *mpcv2+meu+msc+mpv+£c 
if (modtyp.eq.l ) then 

call apdavetpa, mrov, x, amr, amtyp, rovp, rovmaa, apd, aveapd, 
t aavapd, totime, mnapd, mxapd, count, tm, j, pave, 

& paave, fcef) 

else 

minapd - mnapd 
maxapd - mxspd 
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range - maxspd-minspd 
do 30 i - 0, range 
v ■ float ( i+minspd ) 
if ( smtyp. eq.l) then 

call seroin(syst, rovp(i,j), mrov, x, v, pa, smr(i,j), 0 ) 

« * S 6 

rovp(i,j) - (pa+x*0. 08/1000. 0 *(mrov*v) )/( 1 . 0 -x* 

* (0. 088/1000. 0*(smr(l, 1)*v) ) ) 

endif 

rovmas ( i , j ) - rovp( i , j ) * 1 . l* B mr ( i, j )+mrov 
30 continue 
endif 
c 
c 


c ** Output Values ** 

c 

call rovout(pa, j, smr, minspd, maxspd, x, »pcvl, «pcv 2 , meu, 

4 msc, na, systyp, mpv, rovp, rovmas, modtyp, avespd, 

4 savspd, mnspd, mxspd, count, tm, spd) 

print*, 'Do you want to change any of the values (yes or no)?' 
read(5,35)rerun 
35 format(a) 

if ( rerun. eq. 'yes' ) then 
j - j + 1 

smr ( 1 , j ) - smr ( 1 , j- 1 ) 

call newval ( pa , smr(l,j), mnspd, mxspd, x, mpcvl, mpcv 2 , meu, 

4 ®sc, na, modtyp, systyp, smtyp, auxtyp) 

90 to 15 
endif 
stop 
end 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


subroutine rovdatfpa, smr, mnspd, mxspd, x, mpcvl, mpcv 2 , 

4 meu, msc, na, modtyp, systyp, smtyp, auxtyp, 

4 syst, spd, tm, count) 


real pa, smr, x, mpcvl, mpcv 2 , meu, msc, spd(0:50), tm(0:50), 

* mnspd, mxspd ' 

integer na, cartyp, systyp, modtype, smtyp, auxtyp, syst, count 


print*, 'Enter number of cars in the rover excluding the ' 
t ' reactor car . ' 

read( 5 , * )numcar 
print*,' ' 
print*,' ' 

do 1010 i - 1 , numcar 
1000 print*,' ' 
print*,' ' 
printl005,i 

1005 format( ' Enter the type of car number ',il) 
print*, 'Primary Control Vehicle (PCV) - 1' 
print*, 'Experimental Unit (EU) - 2' 

print*, 'Storage and Supply Car (SC) - 3' 
read( 5, * ) cartyp 
if ( cartyp. eq. 1 ) then 
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1010 


1020 


1021 


the mass of the EU (in kg).' 


,3) then 


the mass of the SC (in kg) 


print* , ' ' 

pr int* , ' Enter the mass of the PCV independent of the' 
print*, 'number of Astronauts (in kg).' 
read( 5, * )pcvl 
mpcvl - mpcvl+pcvl 
print* , ' ' 

pfiSt*! 'Enter the mass of the PCV dependent of the' 
pr int* ,' number of Astronauts (in kg/astro.). 
read( 5 , * )pcv2 
mpcv2 - mpcv2+pcv2 

else 

if ( cartyp.eq. 2 ) then 
print*,' ' 
print* , ' ' 
print* , 'Enter 
read( 5 , * )eu 
meu - meu+eu 
else 

if (cartyp.eq 
print* , ' ' 
print* , ' ' 

print*, ’Enter 
read( 5, * )sc 
msc m msc+sc 
else 

print* , ' ' 
print* , ' ' 

print* ,' ***INVALID ENTRY*** 
go to 1000 
endif 
endif 
endif 
continue 
print*,' ' 
print* , ' ' 
print* , 'Enter 
read( 5, * )na 
print* , ' ' 
print* , ' ' 
print* , ' Enter 
read( 5 , * )pa 
print*,' ' 
print*,' ' 
print*, 'Enter 
read( 5 , * )x 
print* , ' ' 

print* ) ' Enter the method of reactor specific mass determination. 

print*, 'User specified - 0' 

print* ,' Program calculated - 1' 

read( 5 , * ) smtyp 

print* , ' ' 

print* , ' ' 

sprint*?' Enter the reactor specific mass (in kg/kWe).' 
read( 5, * ) smr 
print*,' ' 
print* , ’ ' 

elseif (smtyp. eq.l ) then 

orint*. 'Enter the converter type.' . . 

print*,' 1 - Tree Piston Stirling Engine (He working , 

4 print*,' ** U *2*.l closed Brayton Cycle (He/Xe working fluid).' 


the number of astronauts.' 


the power req. 


'd in addition to mobility (in kWe) 


the traction factor modifier x. 
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print*,' 3 - Open Brayton Cycle (C02 working fluid).' 

print*,' 4 - Thermoelectric conversion (SiGe/GaP).' 

read(5,*)syst 

print*,' ' 

print*,' ' 

if (syst.lt.l.or .syst.gt.4 ) then 
print*,' ' 
print*,' ' 

print*, '***INVALID ENTRY***' 
print*,' ' 
go to 1021 

elseif (smtyp.lt.O.or.smtyp.gt.l ) then 
print*,' ' 
print*,' ' 

print*, '***INVALID ENTRY***' 
print*,' ' 
go to 1020 

else 

smr - 0.0 
endif 
end if 

1040 print*, 'Enter the traversing scenario.' 
print* , 'Constant speed - 0' 
print*, 'Variable speed - 1' 
read( 5,*)modtyp 
print*,' ' 
print*,' ' 

if (modtyp.eq.0) then 

print* , 'Would you like to calculate the mass for a' 

print* ,' Single speed - 0' 

print* , 'Range of speeds - 1.' 

read(5,*)systyp 

print*,' ' 

print*,' ' 

1050 if ( systyp.eg. 1 ) then 

print*, 'Enter the minimum Rover speed (in km/hr).' 
read( 5, * )mnspd 
print*,' ' 
print*,' ' 

print* ,' Enter the maximum Rover speed (in km/hr).' 
read( 5, * )mxspd 
print*,' ' 
print*,' ' 

else 

if ( systyp.eg. 0 ) then 

print*, 'Enter the Rover speed (in km/hr).' 
read ( 5 , * )mnspd 
mxspd - mnspd 
print*,' ' 
print*,' ' 
else 

print*,' ' 
print*,' ' 

print*, '***INVAL1D ENTRY***' 
go to 1050 
endif 
endif 

else 

if (modtyp.eg.l) then 

call sendat (count, tm, spd) 

else 

print*,' ' 
print*,' ' 

print*, '***1NVALID ENTRY***' 
go to 1040 
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end if 
•ndif 

print*, 'Enter ' 

print*, '0 - Photovoltaic auxiliary power system' 

print*, '1 - No PV''s in the auxilairy system’ 

read( 5,*)auxtyp 

return 

end 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


subroutine pvfc(na, mpcvl, mpcv2, x, mpv, meu, msc, £c, fceff) 


real mpvd, mpcvl , mpcv2, mpcv, mfcl, mcells, meu, bsc, fc, 
mfc2, mfc3, mpv 


integer i, j, na 

2000 print*, 'This routine is designed to either accept input from' 

print*,' the emergency scenario from a data file called auxdat.' 

print*, 'If this file does file does not exist then you Bay ' 

print* ,' specify a full prompting for this information enter:' 

print*,' 0 - Senario prompting' 

print*,' 1 - Read from Auxdat' 

read(5,*)ptyp 

if (ptyp.eq.O) then 

call auxpromt (sma, tc, fceff, tad, tan# pie# ve# ta# 
t g, pveff, ta2, ta3, apvs) 

print* , tmd, tmn, ta , tc 

else 

if (ptyp.eq.l) then 

open ( uni t»30 , file- 'auxdat' , status- 'old' ) 
rewind 30 

read( 30, * ) sma, tc, fceff, tmd, tan, pis, ve, ta, g, pveff, 

& ta2, ta3 , apvs 

close(unit-30 ) 

else 

print*, '***INVALID ENTRY***' 
go to 2000 
endif 
endif 

td - 24.62 

mcells • 2.0 

mfcl - ta*pls*na*(0.5+mcells) 

apcv - mpcvl+mpcv2*float(na)+mfcl 

smpv » sma/( g*pvef f ) *1000 . 0 

tarn ■ 12.3 

do 2006 j - 0,3 

tarn - float(j)*4.1 

do 2005 i - 12, 0,-1 

tmn float(i)/1.5 

tc - 12 . 3- tmn 

tan * tmn+tmm 

mpvd - ( na*pls* ( tc+ ( td-tc )/f cef f )-12 . 3*apvs*g*pvef f/1000 .+x* 
t ( 0 . 08/1000 . 0 ) * (mpcv+apvs*sma ) *ve* ( tmd+tmn/f cef f ) ) *aapv/ 

£ ( tc-x* ( 0 . 08/1000 . 0 ) *ve* ( tmd+tmn/f cef f ) *smpv ) 

mpv - Bpvd+apvs*sma 
write(40,*)tc, mpv 
2005 continue 
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2006 continue 

mfc2 - ta2*pls*na* ( 0. 5+mcells ) 
mfe3 - ta3*pls*na*(0. 5+mcells) 

£c - mfcl+mfc2+mfc3 
tatot - (ta+ta2+ta3)*pls*na 
write (35,2010) 

write (35,2020)smpv, fceff, tc, apvs, tmd, tmn 

C mfcl, mfc2, mfc3 

'Auxiliary Power System Input',//) 
’PV specific mass 
’Fuel cell efficiency 
’Deployable PV collection time 
’Stationary PV area 
’Mobility time during the day 

’Mobility time during the night 

’Power needed for emergency LS 


pis, ve, tatot, 


2010 format(25x, 
2020 format(15x, 
15x , 
15x, 
15x, 
15x, 
15x, 
15x, 

#/ » 
15x, 
15x, 
15x , 
15x, 
15x, 

return 
end 


Auxiliary cruising speed 
Total fuel cell reserve 
PCV fuel cell mass 
EU fuel cell mass 
SC fuel cell mass 


f 7 . 2 , ' 

kg/kWe ' ,/, 

£7.2,/ 

» 

£7.2, ' 

hr' ,/, 

£7.2, ' 

■**2 ' ,/, 

f 7 . 2 , ' 

hr' ,/, 

£7.2, ' 

hr * , / , 

£7.2, ' 

kWe/Astro' 

f 7 . 2 , ' 

km/hr' ,/, 

f 7 . 2 , ' 

kWh' ,/, 

11.2,' 

fcg' ,/, 

tl .2 

kg' ,/, 

£7.2, ' 

kg' ,///) 


********************************************************************* 


subroutine zeroin(6yst, rovp, mrov, x, v, pa, smpw, m2) 


real rovp, mrov, x, v, pa, kl, k2, k3, ma, step 
integer syst, m2 

double precision cl, c2, cB, c4, c5, c6, cl, c8, c9, clO, ell, 
& cl2, etal, eta2, eta3 

ma -5.0 
sig - 5.67e-8 
f - 0.8 
tamb - 200. 
alph - 0.2 
eps - 0.85 
step - 100. 

kl - pa+x* ( 0 . 08/1000 . *mrov*v) 

k2 - x*0. 088/1000. *v 

print*, kl,k2 
print*, syst 

if (syst.eq.l) then 
trad - 500. 
tr - 2.7 
smc - 6.0 

cl - 0.0543767273338744 

c2 - -0.3981779728928814 
c3 - 0.9497880537528545 

c4 - -0.7477923527039820 
c5 - -0.7487642015184974 
c6 - 5.5454480019398034 

cl —13.4344707466661930 
c8 - 10.7927102688699961 
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c9 - 2.6976636228013784 

clO —20.2721544168889523 
ell - 50.0956199839711189 
cl2 —40.5948684886097908 
c 

•tal - (cl*(tr**3)+c2*(tr**2)+c3*tr)*sac**2 
• ta2 - (c4*smc**2)+(c5*(tr**3)+c6*(tr**2)+c7*tr)*sac 
eta3 ■ (c8*sac)+(c9*tr**3)+(cl0*tr**2)+(cll*tr)+cl2 
ata - etal+eta2+eta3 
•ta ■ ( 1 .-1 ./tr ) *eta 
c 

anc - sac*eta 
elseif (syst.eq.2) then 
•ta - 0.361 
trad - 395. 
sac - 15.75 
elseif ( syst .eq. 3 ) then 
•ta - 0.43 
trad - 1000.0 
sac - 50.0 

•ls«i£( syst.eq. 4 ) then 
•ta - 0.053 
sac - 0.609 
trad - 800. 
endif 

print*, eta, sac 

k3 (l-«ta)*aa/(sig*f*(eps*trad**4-alph*taab**4) )*1000. 
a -4-sac 

print* , k3 
c 

•rr - 0.01 
rl - 10. 
lchk - 1 
nchk - 1 
fr - «rr+l. 
c 

if (m2.eq.l )then 

sapw • 66.54*rovp**(-0.7189)+1494.*rovp**(-0.6203)+k3 
go to 7075 

•lse 

7000 sapw - 66.54*rl**(-0.7189)+1494.*rl**(-0.6203)+k3 
if ( ( k2*sapw/eta) .ge . 1 . ) go to 7050 
fr - kl/( 1 .-k2*sapw/eta)-rl 
print*, rl, fr, step, sapw 
nchk - abs(fr)/fr 
print*, nchk, lchk 
if (nchk. n*. lchk) step - -step/2. 

7050 if (abs(fr) .gt.err ) then 
rl - rl-4-step 
lchk - nchk 
go to 7000 
•naif 
c 

rovp - rl 
•ndif 
c 

7075 return 
•nd 


c 

c 

subroutine rovout (pa, j, sar, ainspd, aaxspd, x, apcvl, apcv2, 
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& 

& 

t 


■«u, msc, na, systyp, mpv, rovp, rovmas, 
modtyp, avespd, savspd, mnspd, mxspd, count, 
tm, speed) 


c 

c 

c 

c 

real pa, smr (0:50,5), npcvl, mpcv2, meu, nsc, mpv, 
a rovp( 0:50,5) , rovmas ( 0 : 50 , 5 ) , avespd, savspd, 
a mnspd, mxspd, tm(0:50), speed(0:50) 
integer j, i, na, systyp, maxspd, minspd, spd, modtyp, count 


write(35,3000) 

3000 format( 22x, 'Hass and Power Requirements for Nuclear',/, 

a 28x, 'Reactor Powered Hars Rover ' ,//,36x, 'Input Values',/) 

if ( modtyp. eq.0) then 

write (35, 3010) pa, minspd, maxspd, x, mpcvl, mpcv2, 
a meu, msc, mpv, na 

3010 format( 15x, 'Supplemental power above mobility - ',f8.2,' kWe' 


a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

3020 

a 

a 


3030 


a 

3031 

3032 


» / » 

15x, 'Rover minimum speed 

15x, 'Rover maximum speed 

15x, 'Rover traction factor 

15x, 'Astronaut independent PCV mass 

15x, 'Astronaut dependent mass 


» i8 , ' km/h',/, 

, i8 , ' km/h',/, 

, f8 . 2 ,/, 

, f 8 . 2 , ' kg',/, 

, f8 . 2 , ' kg/astr ' 


»/ 1 

15x,'Hass of experimental unit 
15x,'Mass of the supply car 
15x,'Mass of the Photovoltaics 
15x, 'Number of astronaunts 


,£8.2, * kg',/, 
, f 8 . 2 , ' kg',/, 
, f 8 . 2 , ' kg',/, 
» i8 ,///) 


write( 35, 3020 ) 

format! 31x, 'Output Values' ,/,13x, 'Rover Speed' , 5x, 'Reactor ', 
'Power' ,3x, 'Rover Hass' ,3x, 'Reactor SH' ,/,16x, 'km/h' ,14x, 'kWe' 
1 1 x , 'kg' ,10x, ' kg/kWt ' ) 
do 3030 i - 0, maxspd-minspd 
sd - i+minspd 

write! 35, 3040 ) sd, rovp(i,j), rovmas(i,j), smr(i,j) 
continue 

else 

write! 35 , 3031 )pa, mnspd, mxspd, avespd, savspd, x, mpcvl, 
mpcv2 , meu, msc, mpv, na 

format! 15x, ' Supplemental power above mobility - ',f8.2,' kWe' 


,/, 

15x, 

15x , 

15x, 

15x, 

15x, 

15x, 

15x, 


Rover minimum speed 

Rover maximum speed 

Rover ave speed per travel time 

Rover ave speed per day speed 

Rover traction factor 

Astronaut independent PCV mass 

Astronaut dependent mass 


9 


* 


9 


f 8 . 2 , ' 
f 6 . 2 , ' 
f 8 . 2 , ' 
f 6 . 2 , ' 
f 8 • 2 ,/, 
f 6 . 2 , ' 
f 8 . 2 , ' 


km/h' ,/, 
km/h' ,/, 
km/h' ,/, 
km/h' ,/, 


kg ' ,/, 

kg/astr' 


,/, 
15x , 
15x, 
15x, 
15x, 


Hass of experimental unit 
Hass of the supply car 
Hass of the Photovoltaics 
Number of astronaunts 


, f 8 . 2 , ' kg',/, 
, f 8 . 2 , ' kg',/, 
, f 6 . 2 , ' kg',/, 
,18, /// ) 


write! 35 , 3032 ) 


t 


format! 30x, 'Rover Speed Power Breakdown' ,//,23x , 'Rover Speed', 
2x, 'Travel Time' ,2x, 'Rover Power' ,/,26x, ' (km/h) ' ,9x, 

' (h) ' ,8x, ' (kWe) ' ) 


do 3033 12 - 1, count 

write (35, 3034) speed! 12 ) , tm( i2 ) , rovp( 12, j ) 

3033 continue 

3034 format! 26x, f 5.2, 9x , f 5 . 2 , 6x , f7 . 3 ) 
write! 35, 3035 ) 
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3035 foraat(//,31x, 'Output Values' ,//,13x, 'Rover Speed' ,2x, 'Ave Re', 
& 'actor Power ' ,2x, 'Rover Hass' ,3x, 'Reactor 5M' ,/,16x, 'km/h' ,14x, 
a 'k' , 'We' ,llx, 'kg' ,10x, 'kg/kwt' ) 

write(35,3045)avespd, rovp( count+1 , j ) , rovmas(l,j), smr(l,j) 
endif 

3040 format(16x,i3,10x,f9.3,7x,f8.2,5x,f8.2) 

3045 foraat(16x,f5.2,8x,£9.3,8x,f8.2,4x,f8.2) 
write(35,3050) 

3050 format <///, 80 < '*' ) ) 
return 
end 


subroutine auxpromt( snpv, tc, fceff, tmd, tan, pis, ve, ta, g, 
a pveff, ta2, ta3, apvs) 


print*,' ' 
print*,' ' 

print*, 'Enter the specific aass of the PV cells (in kg/m**2).' 
read( 5, * )sapv 
print*,' ' 
print*,' ' 

print*, 'Enter the time of reserve power in PCV''s fuel cells', 
a ' ( in hours ) . ' 

read( 5, * )ta 
print*,' ' 
print*,' ' 

print*, 'Enter the time of reserve power in EU' r s fuel cells', 
a ' ( in hours ) . ' 

read( 5, * ) ta2 
print*,' ' 
print* , ' ' 

print* ,' Enter the time of reserve power in SC''s fuel cells', 
a ' ( in hours ) . ' 

read( 5 , * ) ta3 
print*,' ' 
print*,' ' 

print* , 'Enter the efficiency of the Fuel Cell cycle.' 
read( 5, * ) fcef f 
print*,' ' 
print*,' ' 

print* , 'Enter the power needed for life support (in kWe/Astro).' 
read( 5, * )pls 
print*,' ' 
print*,' ' 

print* , 'Enter the daily collection tiae for deployable PV''s', 
a ' (in hours) . ' 

read(5,*)tc 
print*,' ' 
print*,' ' 

print* , 'Enter the area of stationary PV"s (in sq. asters)' 

read(5,*)apvs 

print*,' ' 

print*,' ' 

print* , 'Enter the average solar insolation (in kw/a**2).' 

read(5,*)g 

print*,' ' 
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print*,' ' 

print* , 'Enter the efficiency of the PV panels.' 

read(5,*)pveff 

print*,' ' 

print*,' ' 

print*, 'Enter the time of nobility during the day (in hours).' 

read(5,*)tnd 

print*,' ' 

print*,' ' 

print*, 'Enter the tine of nobility during the night (in hours).' 

read(5,*)tnn 

print*,' ' 

print*,' ' 

P r ^ n t* ,' Enter the speed for energency return (kn/hr).' 
read( 5, * )ve 
print*,' ' 
print*,' ' 

return 

end 






c* ************************************** 
c 




c 

c 

c 


subroutine newval (pa, snr, ninspd, maxspd, x, npcvl, npev2, 

* neu, nsc, na, nodtyp, systyp, sntyp, auxtyp) 


real npcvl, mpcv2, neu, nsc, ninspd, naxspd 

4000 vrite(5,4010)pa, snr, ninspd, naxspd, x, npcvl, npcv2, neu, 
* nsc, na 

4010 fornat( lx, 'Which value would you like to change?' 

£ Cw 9 1 1 4 1 __ a . . _ . » 


£ 

£ 


5x, 

5x, 

9 

5x, 

5x, 

5x, 

5x, 

5x, 

, 

5x, 

5x, 

5x, 

5x, 




/ 


read( 5, * ) nval 


1- Supplemental power above nobility* 

2- Specific nass of the reactor 

9 

3- Rover minimum speed 

4- Rover maximum speed 

5- Rover traction factor modifier 

6- Astronaut independent PCV mass 

7- Astronaut dependent mass 

t 

8- Mass of experinental unit 

9- Mass of the supply car 

10- Nunber of astronaunts 

11 - None' ,//) 


A 

,f8.2, 

* f 8 . 2 , 

, f 8 . 2 , 

, f 8 . 2 , 

, f 8 . 2 ,/, 
, f 8 . 2 , 

, f 8 . 2 , 

,f8.2, 

, f 8 . 2 , 
/i8,A 


kWe' ,/, 
kg/kWe ' 

kn/h' ,/, 
kn/h' ,/, 

kg' ,/, 
kg/astr ' 

kg' ,/, 
kg' ,/, 


if (nval.eq.l) then 

print* ,' Enter the new power value (in kWe).' 
read( 5, * )pa 

? o to 4000 

seif (nval.eq.2) then 

print*,'Ent.r the new reactor specific nass (in kg/kWe). 
r9&u( 5 § * )smr 

? o to 4000 

seif (nval.eq.3) then 

print*, 'Enter the new nininum speed (in kn/hr).' 
read( 5, * )ninspd 
go to 4000 

elseif (nval .eq. 4 ) then 

print* , 'Enter the new naxinun speed (in kn/hr).' 
read( 5, * ) naxspd 
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oo to 4000 

elseif (nval.eq.5) then 

print*, 'Enter the new traction factor modifier.' 
read( 5, * )x 


? o to 4000 

aelf (nval.eq.6) then 
print*, 'Enter the new astronaut 
print*, 'Primary Control vehicle 
read( 5, * )mpcvl 

? o to 4000 

seif (nval.eq.7) then 
print*, 'Enter the new astronaut 
print* ,' Primary Control Vehicle 
read( 5, * )mpcv2 


dependent mass of the' 
(in kg/Astro) . ' 


independent mass of the 
(in kg/Astro) . ' 


? o to 4000 

seif (nval.eq.8) then 

print*, 'Enter the new mass of the Experimental Unit (in kg).' 
read(5,*)meu 


? o to 4000 

seif (nval.eq.9) then 

print* , 'Enter the new mass of the Supply Car (in kg).' 
read(5,*)msc 

? o to 4000 

seif (nval.eq.10) then 
print* , 'Enter the new number of Astronauts.' 
read(5,*)na 


? o to 4000 

seif (nval.ne.il) then 
print*, '***INVALID ENTRY***' 
go to 4000 
endif 


return 

end 

c 

c * ************** ******************************************************* 
C* ********************************************************************* 

c 

subroutine sendat( count, tm, spd) 
c 

real spd(0:50), tm(0:50), tmchk 
integer count 
c 

tmchk - 24.6 
count ■ 1 

print*, 'The length of the Martian day is 24.6 hours. Therefore' 
t , ' , the sum of the ' 

print*, 'time intervals should equal 24.6 hours.' 
c 

6000 print*,' ' 

print*, 'Enter the time interval over which the speed is consent. ' 

read( 5 , * ) tm( count ) 

if ( ( tmchk-tm( count ) ) . It. 0) then 

print* , 'current time interval increases the sum past 24.6 hours' 
print*, 'Tine interval must be less than or equal to ', tmchk 
go to 6000 
endif 

6010 print*, 'Enter the rover speed over the time interval (-> 0km/h).' 
read ( 5 , * ) spd ( count ) 
if ( spd ( count ) .It. 0) go to 6010 
tmchk ■ tmchk-tm( count) 
if ( tmchk. gt.0. 001 ) then 

print*, 'There are still ', tmchk, 'hrs left in the day.' 
count - count+1 
go to 6000 
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or > o nnnn o o o o oo 


endif 


return 

•nd 


subroutine spdave ( pa , mrov , x , smr , smtyp , r ovp , rovmas , spd , avespd , 

* savspd, totime, minspd, maxspd, count, tm, j ,pave, 

& spave,fc) 


real avespd , totime r spd ( 0 : 50 ) , rovp( 0 : 50 , 5 ) ,pa , x , v , 

a rovmas ( 0 : 50 , 5 ) , mrov, smr (0:50,5) , savspd, rp, tm( 0 j 50 ) , 

a maxspd,minspd 

integer i,j, count 

avespd - 0.0 
totime - 0.0 
minspd - spd(l) 
maxspd - spd(l) 
rp - 0.0 


do 8000, i • 1, count 

avespd - avespd+spd( i ) *tm{ i ) 
if ( spd( i ) .ne . 0 . ) totime ■ totime+tm(i) 
print*, totime, tm(i) 
if (spd( i ) .gt. maxspd) then 
maxspd - spd(i) 

elseif (spd(i). It. minspd) then 
minspd - spd(i) 
endif 

6000 continue 
c 

savspd - avespd/24.6 
avespd - avespd/totime 
c 

do 8050 i - 1, 2 
if (i.eg.l) then 
v - avespd 
else 

v ■ savspd 
endif 

if ( smtyp. eg. 1) then 

els* 11 * eroin(syBt ' rov P (i 'i * ,mrov,x,v,pa,smr(i, j ) ,0) 

rovp(i, j ) - (pa+x* ( 0 . 08/1000 . ) *(mrov*v) )/(l . 0-x* 
a ((0. 088/1000. )*(smr(l,j)*v))) 

endif 

if(i.eg.l) then 
pave - rovp(i,j) 

else 

psave ■ rovp(i,j) 
endif 

8050 continue 

print*, pave, psave 

c 

c 

rp - 0.0 

do 8100, i ■ 1, count 
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8100 


c 


v - epd(i) 

If (sntyp.eq.l) then 

call zeroin(ay*t, rovp(i , j ) ,*rov,x,v,pa,a«r ( i , J ) ,0) 

#1 rovp(i,1) - (pa+x*(0. 08/1000. )M»rov*v) )/(1.0-x* 
r ° V P' ^ ( ( 0 . 088/1000. )*(snr (1, j )*v) ) ) 

if Uovp( i , j ) . gt .paave ) covp( i , j )-( rovp( i , j )-p«ave )/fc+paave 
rp ■ ( rovp( i , j )-pa) *tm( i )+rp 
print*, rovp(i, j) 
continue 


rovp( count+1 


,j) - rp/24.6+pa 

if T amtyp.eq. 1 ) then 

call zeroin( sy*t , rovp( count+1 , j ) ,»rov,x, v ,pa , anr ( 1 , j ) , 1 ) 
rovmas(l,j) - rovp( count+1 , j )*l.l*»mr(l, j )+ar° v 


return 

end 
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APPENDIX C 

VALUES OF PHYSICAL CONSTANT 
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•Prom L. T. Templin. editor. Reactor Phytict Conel.nl*. ANL5800, 2nd ed.. 1963; baaed on G. W. Grodttoin. National Bureau of Standard, circular 583. 1957. 



APPENDIX D 

Sample TWODANT Input File 
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2 0 0 

CORE MODEL FOR MANNED MARS ROVER. 12-1-89 
ADAPTED FROM INPUT CREATED BY JOHN MCGHEE. 


block 1 - controls and dimensions ************ 


IGEOM- 

NGROUP- 

ISN- 

NISO- 

MT- 

NZONE- 

IH- 

IT- 

JM- 

JT- 

MAXLCM- 

MAXSCM- 

NOGEOD- 

NOMIX- 

NOASG- 

NOMACR- 

NOSLNP- 

NOSOLV- 

NOEDIT- 

T 


7 

42 

8 

60 

11 

27 

10 

27 

14 

32 

281500 

18500 

1 


/ r/z geometry 
/ 4 of rieutron groups 
/ Sn order 

/ # of isotopes on xslib 

/ # of materials 

/ I of model zones 

/ I of r coarse mesh intervals 

/ # of r fine mesh intervals 
/ I of z coarse mesh intervals 
/ # of z fine mesh intervals 
/ maximum laroe core memory 
/ maximum small core memory 
/ 1 » suppress geometry input 

/ 1 - suppress mixino input 

/ 1 - suppress material zone assignment input 

/ 1 « suppress genration of macroscopic xs 

/ 1 - suppress generation of solver input 
/ 1 - suppress solver module 

/l - suppress edit module 


block 2 - geometry details 


/ r coarse mesh point coordinates (cm); 

XMESH- 0 3 8 13 16 17 19 23 27 30 230 
/ z coarse mesh point coordinates (cm); 

YMESH- 0 50 53 60 65 73 76 110 113 130 132 133 147 150 200; 
/ fine intervals per r coarse interval: 

XINTS- 122211222 12 
/ fine intervals per z coarse interval: 

TINTS- 52221242211125 
/ zone number for each coarse mesh interval: 

ZONES- 14 14 14 14 14 14 14 14 14 14; 

27 27 27 27 27 27 27 27 27 14; 

1 26 26 26 26 26 26 26 27 14; 

1 2 4 2 3 0 0 0 27 14; 


1 2 4 2 3 

6 2 5 2 3 

6 7 8 7 7 

6 9 10 11 12 

6 18 22 20 12 

17 18 22 20 12 

17 21 24 21 21 

17 23 24 25 25 

17 23 23 25 0 


0 0 27 14; 

0 d 27 14; 
15 13 27 14; 
15 13 27 14; 
15 li 27 14; 


d 27 14 
0 27 14 
0 27 14 
0 27 14 


27 27 27 27 27 27 27 27 27 14; 
14 14 14 14 14 14 14 14 14 14; 


LXB- BXSLXB 

mXxord- 1 

IHT- 4 

LNG- 30 

ITITL- 1 

SAVBXS- 0 

KKIKRD- 0 

/ names for each isotope: 
NAMES- HI D2 T3 


block 3 - riuclear data *************** 

/ cross section source 
/ highest legendre order 
/ total xs row I 
/ lait neutron group 
/ title cards included 
/ 1 - save binary xs file 
/ 1 - no error checking 
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Bll 

C12 

N14 

016 

ri9 

CL 

K 

CA 

TI 

V 

NI 

CU 

2R 

HB93 

MO 

BA138 

TA181 

W182 

W183 

W184 

PB 

TH232 

PA233 

U233 

U234 


NA23 MG AL27 SI 
CR HN55 rE C059 
RH103 AG107 AG109 CD 


ii r 4 j i ruoo ruts* ru^«u rUifll rUZ4Z 

IN EBOUND ARRAY THE LAST NEUTRON GRP LOWER BOUND IS .139-9 MEV 


1.700+1 

1.500+1 

1.350+1 

1*200+1 

1.000+1 

7.790 6.070 

3.680 

2.865 

2.232 

1.738 

1.353 

0.623 0.500 

0.303 

0.184 

0.676-1 

2.480-2 

9.120-3 

3.350-3 1.235-3 

4.540-4 

0.414-6 

1.670-4 

0.152-6 

6.140-5 

2.260-5 

8.320-6 

3.060-6 1.130-6 

10.00 

9.00 

8.00 

7.00 

6.00 

5.50 4.00 

3.00 

2.00 

1.00 

0.50 

0.10 

0.01 

54.8 51 

.9 49.6 

45.0 

40.8 

36.1 

29.4 24.8. 

21.9 19 

.3 17.1 

14.2 

11.1 

8.67 

6.79 4.83 

2.95 1 

.81 1.07 

0.63 

0.39 

0.237 

0.144 0.087 

0.053 0 

for each 

.032 0.020 
edit xs: 

0.012 

0.007 

0.003 

12R300 


S- nss 


******** block 4 - mixing details ********< 
densities (atoms/b-cm) for each material: 


FUEL97 

U238 

.0009564 

U235 

.0309236 

N14 

FUEL73 

U238 

.0086080 

U235 

.0232720 

N14 

LINER 

RE187 

.0412847 

RE185 

.0246653 


PWC11 

NB93 

.0525700 

ZR 

.0005414 

C12 

B4C 

B10 

.0684000 

Bll 

.0076000 

C12 

BEO 

BE9 

.0676600 

016 

.0676600 


BE 

BE9 

.1196000 




BAFFLE 

TA181 

.0540800 




CLNT 

LI 6 

.0027500 

LI 7 

.0343100 


LIH 

HI 

.059 

LI 6 

.0045500 

LI? 

TUNG 

W182 

.0171600 

W18 3 

.0093300 

W184 

ial volume fractions in 

each 

zone : 


■ 1 PWCll 

.190 

; 




2 CLNT 

1.0; 





3 PWCll 

.350 

CLNT .650 

; 



4 PWCll 

.036 

CLNT .800 

; 



5 PWCll 

.038 

CLNT . 800 

BEO 

.134; 


6 PWCll 

.190 

BEO .670 

; 



7 PWCll 

.800 

CLNT .200 





.0318800; 

.0316800; 

.0004112; 

.0190000; 

l 

> 

} 

; 

.0561200; 


W186 .0186700; 


8 PWC11 .678 CLNT 


9 

PWCll 

.129 

CLNT .176 

FUEL97 

.594 LINE& .066; 

10 

PWCll 

.141 

CLNT .141 

FUEL97 

.475 LINER .053 BEO 

11 

PWCll 

.129 

CLNT .176 

FUEL97 

.594 LINER .066; 

12 

PWCll 

.250 

CLNT .650 

BAFFLE 

.100; 

13 

BEO 

.7 

B4C .3; 



14 

LIH 

1.0; 




15 

BEO 

.550 

B4C .450; 



16 

BEO 

.878 

t 



17 

PWCll 

.190 

B4C .670; 



18 

PWCll 

.129 

CLNT .176 

LINER . 

107; 

19 

PWCll 

.141 

CLNT .141 

LINER . 

086 BEO .134; 

20 

PWCll 

.129 

CLNT .176 

LINER . 

107; 

21 

PWCll 

.200 

CLNT .800; 



22 

PWCll 

.678 

CLNT .160 

B4C 

134; 

23 

CLNT 

1.0; 




24 

PWCll 

.038 

CLNT .800 

B4C 

134; 

25 

PWCll 

.250 

CLNT .750; 



26 

BE 

1.0; 
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27 TUNG 1.0; 


T 

/ ************** block 5 — solver details ******************i** 

/ 


IEVT- 

i 

/ eigenvalue search 

ISCT- 

0 

/ scattering legehdre order 

XBL- 

1 

/ left boundary reflective 

IBR- 

0 

/ right boundary vacuum 

IBT- 

0 

/ top boundrary vkcuum 

XBB- 

0 

/ bottom boundary vacuum 


T 

/ 

/ ************ block 6 - edit input details ********************** 

/ 

/ spectrua data 

/ 

PTED- 1 
ZNED- 0 
POINTS- 30127,664 
/ POINTS- 2511,27 
/ POINTS- 251838,864 
/ POINTS- 3, 432, 840 
/ 

/ edit location: lower plenum, 

/ fine intr (r,s): (11,14), 

/ edit points - r+44(z-l): 

/ 

RESDNT- 1 
IGRPED- 2 
ICOLL- 30, 12 
POWER- 0.5 
HEVPER- 200.0 
RSFE- F1.0 
RZFLUX-0 
T 


/ edit by mesh point 
/ no zone edits 


lecore, hecore, rod plenum, upper plenum 
(6,40), (16,40), (11,70), (9,94) 


/ use edit xs resident macro xs 
/ print broad groups and totals 
/ energy group collapsing 
/ normalise output to this value (mw) 
/ mev per fission 
/ responce function multiplier 
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APPENDIX E 

Tungsten Nuclear Cross Section Data 
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l .0"10 ' s .0010 .010 .10 1.0 10.0 100. 1.0-10 3 1.0 x 10 s 2.0-10 

E„ (eV) 


